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Fig.  1.  Acoustooptic  synthetic  aperture  radar  processor  utilizing  discrete  components  6 
[1].  The  top  view  shows  the  range  focus  due  to  the  linear  fm  SAR  return;  the  side 
view  shows  the  vertical  expansion  in  a  given  range  bin  required  for  time  integration 
on  the  CCD  array  detector. 

Fig.  2.  Schematic  diagram  of  the  Integrated  Optical  Synthetic  Aperture  Radar  Pro-  10 
cessor,  as  described  in  Sect.  4  of  the  text. 

Fig.  3.  Waveguide  outcoupling  from  a  rib  or  channel  waveguide  by  uniform  periodic  15 
modulation  of  the  waveguide  surface. 

Fig.  4.  Waveguide  outcoupling  from  a  rib  or  channel  waveguide  by  nonuniform  periodic  16 
modulation  of  the  waveguide  surface  to  compensate  for  outcoupling  loss  in  the 
direction  of  propagation. 
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The  principal  objective  of  this  research  program  was  initiation  of  the  development  of 
a  compact,  rugged,  low-cost,  low-power  optical  synthetic  aperture  radar  processor  for 
real-time  image  formation  aboard  airborne  and  spaceborne  platforms. 

The  extremely  high  computational  burden  associated  with  SAR  imaging  has  thus  far 
limited  the  usefulness  of  the  technique  to  non-real-time  applications  from  fixed  platform 
bases.  Existing  real-time  digital  and  optical  techniques  are  in  general  too  bulky,  power 
consumptive,  and  expensive  for  broad  applicability,  particularly  in  space-confined  mo¬ 
bile  environments.  The  integrated  optics  approach  described  herein  provides  a  powerful 
solution  with  which  to  circumvent  these  limitations. 

The  key  to  this  novel  optical  processing  concept  is  the  use  of  planar  optical  wave¬ 
guide  technology  (which  allows  the  integration  of  optical  components  in  a  single  sub¬ 
strate)  in  conjunction  with  a  time-and-space  integrating  architecture  (which  allows  two- 
dimensional  processing  to  be  performed  in  a  quasi-planar  optical  structure).  An  addi¬ 
tional  novel  feature  is  the  utilization  of  partial  waveguide  confinement  to  allow  uniform 
light  emission  vertically  out  of  the  plane  of  the  integrated  optical  substrate  and  onto  a 
two-dimensional  time-integrating  detector.  ^ 

The  unique  features  of  this  approach  are:  / 

1.  The  implementation  of  two-dimensional  SAR  processing  utilizing  both  integrated 
optics  and  time-and-space  integration 

2.  Vertical  integration  of  an  optical  waveguide  with  a  CCD  image  sensor 

3.  The  use  of  in-plane  parallel  channel  or  rib  waveguiding  to  retain  lateral  image 
integrity 

4.  Surface  acoustic  wave  generation  in  optical  waveguide  materials  with  slow  acoustic 
velocity  to  achieve  large  effective  radar  apertures 

5.  Selective  partial  waveguide  confinement  of  channelized  light  in  an  integrated  optics 
chip 
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The  novel  integrated  optics  architecture  described  in  this  report  is  potentially  pro¬ 
grammable  to  allow  additional  types  of  computational  algorithms  to  be  implemented. 
These  include  the  compensation  of  SAR  images  for  changes  in  platform  variables  and 
geometric  distortions,  two-dimensional  correlation,  convolution,  and  Fourier  transform 
operations,  and  real  time  computation  of  ambiguity  functions. 
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1.  INTRODUCTION 

Synthetic  Aperture  Radar  (SAR)  is  a  very  powerful  imaging  technique  that  allows  the 
formation  of  high  resolution  images  using  radar-frequency  probes.  In  a  conventional  real 
aperture  imaging  system,  the  resolution  is  inversely  proportional  to  the  aperture  of  the 
imaging  system.  Thus,  improved  resolution  can  only  be  obtained  at  the  expense  of  a  large 
radar  antenna,  which  is  in  general  too  bulky  to  carry  aboard  an  airborne  or  spaceborne 
platform.  In  many  cases,  the  desired  resolution  would  require  antennas  of  truly  massive 
size,  which  are  impractical  for  any  platform.  In  a  SAR  system,  on  the  other  hand,  the 
resolution  is  approximately  equal  to  the  antenna  aperture  in  the  azimuth  direction  and 
it  is  entirely  independent  of  the  antenna  aperture  in  the  range  dimension. 

The  signal  detected  by  a  SAR  system,  however,  is  not  a  focused  image;  extensive 
processing  is  required  to  form  the  image  from  the  received  signal.  Film  based  optical 
processing  systems  have  been  very  effectively  used  in  the  past  twenty-five  years  to  per¬ 
form  the  required  SAR  processing.  More  recently,  digital  computers  have  been  also  used 
extensively  in  SAR  image  formation.  However,  film  based  optical  processors  and  high 
speed  minicomputers  are  not  suitable  for  real  time  image  formation  aboard  the  radar 
platform  because  of  computational  speed  limitations,  and  also  because  of  the  large  size, 
weight  and  power  consumption  of  such  systems.  Recently,  optical  techniques  have  been 
demonstrated  at  both  Caltech  and  USC  that  allow  real  time  SAR  image  formation  with 
relatively  compact  and  power  efficient  optical  systems.  These  systems  were  implemented 
with  discrete  components  [1,2,3;  see  also  Appendices  A,  B,  C).  In  one  approach  [1,2], 
an  acoustooptic  cell  was  employed  to  input  each  radar  return  into  a  time-and-space 
integrating  optical  architecture  comprised  of  severed  lenses,  a  CCD  area  array  detec¬ 
tor,  emd  a  proximity-coupled  mask.  In  this  approach,  the  range  compression  derives 
from  self-focusing  of  the  chirped  radar  pulse,  while  azimuth  compression  is  achieved  by 
mask  multiplication  and  sequential  time-integration  on  the  CCD  array,  operated  in  the 
shift-and-add  mode.  This  architecture  is  explained  in  more  detail  in  Section  2  and  in 
Appendices  A  and  B.  In  the  second  approach  [3],  the  incoming  radar  return  at  a  given 
azimuth  coordinate  is  written  onto  one  column  of  a  two-dimensional  real-time  spatial 
light  modulator  such  as  the  PROM  (Pockels  Readout  Optical  Modulator)  [4],  Sequential 
radar  return  pulses  are  written  onto  adjacent  columns.  This  input  function  can  be  per¬ 
formed  by  a  serial-to-parallel  converter  such  as  the  Linear  Array  Total  Internal  Reflection 
Spatial  Light  Modulator  [5].  Range  and  azimuth  compression  are  then  both  achieved  in 


parallel  by  means  of  an  anamorphic  imaging  system,  which  relies  on  the  self-focusing 
properties  of  the  radar  return  Fresnel  zone  plate  in  both  dimensions.  Further  details  are 
provided  in  Appendix  C. 

A  further  reduction  in  the  size,  weight,  power  and  cost  requirements  of  a  real-time 
SAR  processor  can  be  realized  if  these  optical  components  could  be  integrated  onto  a 
single  optical  waveguide  substrate.  The  value  of  such  an  achievement  can  be  immediately 
realized  by  comparison  of  several  pertinent  characteristics  of  existing  real-time  SAR 
processors.  For  example,  a  current  digital  implementation  comprises  the  equivalent  of 
two  vertical  six-foot  equipment  racks,  weighs  about  1500  pounds,  consumes  from  one 
to  five  kilowatts  of  power,  and  costs  several  hundred  thousand  dollars  [6].  The  optical 
implementations  described  above  are  of  order  one-half  cubic  meter  in  size,  weigh  about 
twenty  pounds,  consume  less  than  five  watts  of  power,  and  cost  about  ten  thousand 
dollars.  The  integrated  optical  approach  described  herein  will  be  of  order  four  cubic 
inches  in  size,  weigh  about  one-half  pound,  consume  less  than  200  milliwatts  of  power, 
and  cost  no  more  than  a  few  hundred  dollars  in  quantity. 

The  novel  aspect  of  the  proposed  research  is  the  computation  of  two-dimensional 
operations  using  integrated  optics.  The  formation  of  a  SAR  image  is  inherently  a  two- 
dimensional  operation,  in  that  both  azimuth  and  range  information  must  be  extracted 
from  the  microwave  radar  returns.  Integrated  optics,  on  the  other  hand,  is  a  planar 
technology  usually  thought  suited  for  only  one-dimensional  processing  functions  (and  to 
which  it  has  been  recently  applied  [7,8,9].  In  this  report,  we  describe  initial  feasibility 
studies  of  the  implementation  of  two-dimensional  processing  with  an  integrated  optics 
structure,  using  the  time-and-space  integrating  architecture  [10].  With  this  method, 
two-dimensional  data  is  encoded  in  time  and  one  spatial  coordinate,  and  the  data  is 
processed  by  a  combination  of  spatial  and  temporal  integration  within  the  optical  system 
to  achieve  the  requisite  range  and  azimuth  compression  functions.  Although  only  one 
spatial  coordinate  is  utilized,  it  is  possible  to  accomplish  extremely  complicated  two- 
dimensional  signal  processing  functions  such  as  SAR  image  formation  within  the  confines 
of  a  planar  integrated  optics  structure  by  means  of  a  novel  architecture  that  involves  the 
vertical  extraction  of  image  intensity  from  an  array  of  parallel  waveguides.  Temporal 
integration  is  provided  by  a  proximity-coupled  area  array  CCD  detector.  This  concept 
is  described  more  fully  in  Section  4,  Summary  of  Research  Program. 


2.  BACKGROUND  OF  THE  RESEARCH  PROGRAM 


A  schematic  diagram  of  the  aforementioned  discrete  acoustooptic  SAR  processor  [1,2] 
is  shown  in  Fig.  1.  The  theory  of  operation  of  this  processor  is  described  in  detail 
in  Appendix  A,  and  a  recent  experimental  demonstration  is  presented  in  Appendix  B. 
Briefly,  the  input  to  the  processor  is  the  radar  signal  return,  which  is  applied  to  the 
acoustooptic  device  (AOD)  in  Fig.  1.  The  light  diffracted  by  the  AOD  is  expanded 
uniformly  in  the  vertical  direction  and  it  self-focuses  horizontally  due  to  the  linear  FM 
pulse  shape  of  the  radar  signal.  The  self-focusing  of  the  diffracted  light  produces  a 
focused  radar  image  in  the  range  dimension  only.  The  range  focused  image  is  multiplied 
by  a  fixed  two-dimensional  mask,  on  which  the  azimuth  phase  history  is  recorded,  and 
the  product  is  detected  on  a  two-dimensional  CCD  detector  array.  The  photogenerated 
charge  on  the  CCD  is  continuously  transferred  step-wise  in  the  vertical  direction  during 
the  exposure  of  the  device  to  light.  The  relative  motion  between  the  signal  on  the  CCD 
and  the  stationary  mask  results  in  the  required  correlation  in  the  azimuth  direction. 
Thus  the  output  video  signal  from  the  CCD  is  the  two-dimensional,  focused  radar  image. 
In  this  architecture  the  simultaneous  use  of  time  and  space  integration  allows  the  two- 
dimensional  image  forming  operation  to  be  implemented  with  a  one-dimensional  AOD 
as  the  input  transducer.  This  implementation  utilizes  two  transverse  spatial  coordinates; 
the  second  dimension  is  needed  to  uniformly  expand  the  light  in  order  to  address  the 
two-dimensional  mask  and  CCD.  In  an  integrated  optics  (10)  implementation  of  such 
a  processor,  a  method  must  therefore  be  found  for  addressing  the  mask  and  CCD  in  a 
planar  geometry.  This  method  is  described  in  Section  4. 

The  term  ’’integrated  optics”  refers  to  the  technology  with  which  optical  and  electro- 
optical  components  can  be  integrated  on  a  single  substrate.  The  components  that  we 
are  interested  in  integrating  are  the  light  source,  detector,  lenses,  waveguides,  and  mod¬ 
ulators.  Therefore,  in  selecting  a  substrate  material,  consideration  must  be  given  to  the 
technology  that  is  available  for  fabricating  each  of  the  individual  components  with  this 
material.  At  present,  no  one  material  exists  within  which  all  the  required  components  can 
be  fabricated.  Consequently,  fully  integrated  optical  systems  are  not  yet  feasible.  The 
material  most  commonly  used  thus  far  is  LiNb03  ,  in  which  optical  waveguides,  lenses, 
electrooptic  and  acoustooptic  (SAW)  modulators  can  be  monolithically  fabricated.  The 
laser  diode  source  and  detector  are  typically  butt-coupled  at  the  two  ends  of  the  inte¬ 
grated  optics  chip.  Materials  such  as  GaAs  show  promise  for  full  integration  in  future 
years. 

To  date,  integrated  optical  (or  at  least  hybrid-integrated  guided  wave  optical)  spec¬ 
trum  analyzers,  one-dimensional  correlators,  and  simple  vector-matrix  multipliers  have 
been  demonstrated.  One  of  the  frequently  stated  limitations  of  integrated  optics  technol- 
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Fig.  1  Acoustooptic  synthetic  aperture  radar  processor 
utilizing  discrete  components  [1].  The  top 
view  shows  the  range  focus  due  to  the  linear  fm 
SAR  return;  the  side  view  shows  the  vertical 
expansion  in  a  given  range  bin  required  for  time 
integration  on  the  CCD  array  detector. 


3.  STATEMENT  OF  TECHNICAL  OBJECTIVES 


The  full-scale  development  of  an  Integrated  Optical  Synthetic  Aperture  Radar  Proces¬ 
sor  involves  a  wide  range  of  definable  tasks,  including  design,  fabrication,  and  character¬ 
ization  of  a  hybrid-integrated  laser  diode,  a  waveguide  surface  acoustic  wave  transducer, 
several  integrated  optical  collimation  lenses,  a  rib  or  channel  waveguide  array,  SAR  az¬ 
imuth  compression  masks,  a  charge-coupled  device  (CCD)  to  waveguide  interface,  system 
integration,  device  characterization,  demonstration  of  real-time  SAR  image  formation, 
and  the  exploration  of  applications  to  other  significant  problems  amenable  to  solution  by 
optical  processing  algorithms.  In  order  to  accomplish  these  tasks,  a  collaborative  effort 
between  research  groups  at  the  University  of  Southern  California  (Principal  Investigator: 
Professor  Armand  R.  Tanguay,  Jr.)  and  the  California  Institute  of  Technology  (Prin¬ 
cipal  Investigator:  Professor  Demetri  Psaltis)  was  proposed,  and  implemented  by  two 
separate  research  contracts.  In  general,  the  source  hybridization,  surface  acoustic  wave 
device  design  and  fabrication,  and  SAR  algorithm  mask  design  were  designed  to  be  the 
responsibility  of  Caltech,  while  the  integrated  optical  lens  design  and  fabrication,  rib  or 
channel  waveguide  design  and  fabrication,  and  CCD  detector  interface  were  designed  to 
be  the  responsibility  of  USC.  System  integration,  device  characterization  and  analysis, 
and  real  time  SAR  image  formation  tasks  were  proposed  to  be  shared  by  both  insti¬ 
tutions,  as  were  the  investigation  of  other  optical  processing  applications  of  the  novel 
architecture. 

Our  specific  goal  for  this  initial  15  month  program  (covering  the  research  period 
from  9/1/84  through  11/30/85)  was  the  demonstration  of  the  feasibility  of  the  approach, 
and  solution  of  several  key  issues  regarding  its  potential  implementation.  In  Section  4, 
the  specific  research  task  items  that  were  pursued  during  the  course  of  the  research  are 
described. 


4.  SUMMARY  OF  RESEARCH  PROGRAM 

Description  of  the  Processor 

The  proposed  Integrated  Optical  (10)  SAR  processor  is  shown  schematically  in  Fig. 
2.  It  consists  of  a  single  crystal  substrate,  the  surface  of  which  has  been  modified  to 
change  the  index  of  refraction  of  a  thin  layer  near  the  surface  and  thus  form  a  planar 
waveguide.  Light  from  a  semiconductor  laser  (attached  to  the  10  chip)  is  butt-coupled 
at  the  left  and  is  guided  in  the  vertical  direction  along  the  surface  of  the  chip.  In  the 
horizontal  direction  the  light  is  diverging,  and  waveguide- lens  LI  is  used  to  collimate 
it.  A  surface  acoustic  wave  (SAW)  is  launched  along  the  top  surface  of  the  chip  in  a 
direction  approximately  perpendicular  to  the  direction  of  propagation  of  the  light.  The 
SAW  is  modulated  by  the  radar  signal  and  causes  a  portion  of  the  incident  light  to  be 
diffracted.  Lens  L2  focuses  the  diffracted  light  at  plane  PO.  The  one  dimensional  light 
distribution  at  plane  PO  is  the  radar  image  focused  in  the  range  dimension  only.  An 
array  of  rib  waveguides  is  fabricated  on  the  10  chip  after  plane  PO.  The  rib  waveguides 
serve  a  two-fold  purpose.  First,  they  retain  the  integrity  of  the  range  focused  image  after 
plane  PO  (i.e.  they  prevent  the  image  from  going  out  of  focus  as  it  propagates  away  from 
plane  PO).  Second,  they  expand  the  range  focused-image  in  the  plane  of  the  chip  so  that 
the  two-dimensional  mask  and  CCD  that  are  part  of  the  acoustooptic  SAR  processor  can 
be  addressed.  This  is  accomplished  by  making  the  top  surface  of  rib  waveguides  leaky, 
so  that  a  portion  of  the  light  escapes  upwards  as  it  travels  in  each  waveguide.  Several 
methods  for  achieving  this  are  discussed  later.  The  mask  and  CCD  are  then  placed  in 
contact  with  the  surface  of  the  10  chip  so  that  the  escaping  light  will  be  multiplied  by  the 
transmittance  of  the  mask  and  detected  by  the  CCD.  The  radar  image  is  focused  in  the 
second  (azimuth)  direction  by  transferring  and  integrating  the  photogenerated  charge  on 
the  CCD  precisely  as  in  the  the  acoustooptic  SAR  processor  described  in  Appendices  A 
and  B. 

A.  Surface  Acoustic  Wave  Device  Fabrication 

SAW  devices  have  been  fabricated  in  a  wide  variety  of  crystals.  Most  of  the  SAW /IO 
devices  are  fabricated  with  LiNbOa  at  the  current  time,  primarily  because  optical  wave¬ 
guides  and  other  optical  components  can  also  be  fabricated  easily  in  LiNbOa  .  However, 
LiNbOa  is  not  optimum  for  the  SAR  application  because  the  acoustic  velocity  is  too 
high,  and  this  limits  the  allowable  acoustic  delay  (and  hence  the  range  delay  in  the  SAR 
processor)  to  approximately  1  microsecond.  In  the  SAR  application,  an  acoustic  delay 
time  of  at  least  50  microseconds  is  required  in  order  to  generate  sufficient  range  coverage. 
We  investigated  the  properties  of  a  range  of  single  crystal  substrate  materials  with  low 
acoustic  surface  wave  velocities,  to  determine  those  in  which  both  SAW  devices  and  opti- 


Ui  WUlWWlWWUTWlPffVJWyvmiV^wnf^wm^^ti^p^y^ 


10 


fO 

u 


CL  CO 

O  *3 


*o  * 

CJ  S- 
4-J  o  • 
TO  cn  +j 

£-  </i  x 
cn  <u  a> 

Q)  UP 
4->  O 
C  S-  a) 
►— *  Q_  _C 


uiolj 


Sglfi 

CO<  £  i.CO 


<U  1- 


^  CO 


So1 


+j  -a  o 

*  •'« 

ra 

•  *  • 

*  *  * 

<4-  oc  *r 

•  V 

o 

Q)  • 

E  j.  +j 

k!kA 

1 

<0  3  0 

r-  * 

S-  4->  CD 

cn  s_  go 

<o  <u 

■c'J 

■<-  a.  c 

.\V 

x>  C  ■<- 

.*  V 

0  0X3 

■A'-' 

•<-  •<-  a> 

AV 

+->+->  JO 

<g  a)  •<- 

E  -C  s_ 

<u  -o  o 

,4 

£  C  «l 

H 

o  >>  cu 

GO  GO  TO 

V? 
■  *  * 

CM 

> 

:-:M 


•:%i 


yn 


i 

/ 


$ 


£s 


.  a'..'  v"-.'' 


cal  waveguides  can  be  fabricated.  Although  tellurium  dioxide  (TeOj)  has  a  significantly 
lower  surface  acoustic  wave  velocity  (shear  wave)  than  that  of  LiNb03,  current  technol¬ 
ogy  does  not  seemingly  permit  straightforward  and  reproducible  fabrication  of  large  area 
waveguides  in  this  material.  This  appears  to  be  primarily  due  to  a  lack  of  investigative 
effort  moreso  than  to  any  fundamental  reason. 

As  a  direct  result  of  this  survey,  it  was  determined  that  LiNbC>3  is  a  more  attractive 
material  for  purposes  of  concept  demonstration.  As  such,  we  decided  to  fabricate  all 
proof  of  concept  device  structures  and  substructures  on  LiNbC>3  substrates. 

B.  Integrated  Optical  Lens  Fabrication 

A  second  key  component  of  the  Integrated  Optical  SAR  Processor  involves  the  fabri¬ 
cation  of  a  collimating  lens  for  the  laser  diode,  and  a  focusing  lens  for  the  acoustooptic 
device  (to  significantly  shorten  the  focal  lengths  that  characterize  typical  linear  FM  chirp 
waveforms  in  a  SAW  transducer).  A  large  number  of  investigators  have  examined  the  lens 
fabrication  issue,  and  a  range  of  implementable  solutions  exists,  including  both  geodesic 
and  Luneberg  lenses.  Since  this  component  has  been  thoroughly  investigated  by  others, 
we  proposed  to  postpone  examination  of  the  specific  issues  involved  in  such  a  lens  de¬ 
sign  until  all  other  critical  issues  that  bear  on  proof-of-concept  have  been  satisfactorily 
demonstrated.  During  the  initial  program  period,  we  did  however  search  the  litera¬ 
ture  for  previous  results  that  bear  directly  on  our  task,  examined  questions  related  to 
lens  f-number  that  derive  from  appropriate  free  propagation-to-confined  waveguide  mode 
propagation,  and  incorporated  the  results  of  these  studies  into  design  of  the  optimum 
acoustooptic  cell  and  parallel  rib  waveguide  structure. 

During  the  course  of  the  literature  survey,  we  became  aware  of  an  elegant  and  poten¬ 
tially  profound  advance  in  the  fabrication  of  integrated  lenses  and  lens  arrays  in  LiNb03 
by  the  titanium  indiffused-  proton  exchange  (TIPE)  process  [11].  Nearly  diffraction  lim¬ 
ited  performance  was  obtained  in  several  experimental  demonstrations  of  lithographically 
defined  microlenses  and  microlens  arrays,  including  lenses  of  large  aperture  [11].  These 
results  are  of  immediate  significance  to  our  current  research,  as  they  point  to  greatly  sim¬ 
plified  but  high  performance  large  aperture  lenses  as  required  by  the  IOSAR  processor 
geometry.  In  particular,  the  extraordinary  large  refractive  index  changes  induced  by  the 
TIPE  process  imply  much  lower  curvature  lens  surfaces  for  the  same  f-number,  which  in 
turn  reduces  the  aberrations  as  well  as  the  overall  size  of  the  processor. 

C.  Rib/Channel  Waveguide  Fabrication 

As  shown  in  Fig.  2,  an  essential  feature  of  the  integrated  optical  SAR  image  formation 
concept  is  the  fabrication  of  a  parallel  array  of  either  rib  or  channel  waveguides  that  si¬ 
multaneously  preserve  the  range  information  while  allowing  the  azimuth  correlation  to  be 
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p>  :  formed  by  the  shift-and-add  CCD  detector.  A  number  of  key  issues  pertain  to  the  suc¬ 
cessful  implementation  of  such  a  novel  structure  (unique  primarily  in  that  the  individual 
waveguides  are  to  be  designed  to  couple  a  given  amount  of  light  per  unit  length  out  of  the 
top  of  the  guide,  through  the  proximity-coupled  mask  and  onto  the  CCD  array  detector). 
These  issues  include  the  minimization  of  cross-talk  between  individual  ribs  or  channels, 
the  choice  of  either  ribs  or  channels  for  array  implementation,  the  coupling  efficiency 
from  the  mesa  containing  the  integrated  lens  into  each  individual  rib/channel  (each  of 
which  represents  an  independent  range  bin),  the  minimization  of  sidewall  irregularities 
that  can  lead  to  scattered  light  loss  (and  a  consequent  reduction  in  both  outcoupling 
uniformity  and  overall  signal-to-noise  ratio),  and  the  uniform  fabrication  of  large  area 
arrays  (of  order  lxl  cm2)  with  very  large  numbers  of  independent  ribs  or  channels  (of 
order  1000).  Further,  the  fabrication  process  chosen  must  be  compatible  with  a  method 
of  outcoupling  from  the  array  that  can  be  produced  with  minimal  degradation  of  the 
array  channel  separation  characteristics.  Methods  of  waveguide  outcoupling  explored 
during  the  contract  period  are  described  in  the  next  Section  (Waveguide  Outcoupling 
Techniques). 

The  first  step  in  the  development  of  a  waveguide  array  is  the  successful  implementa¬ 
tion  of  uniform,  low  loss  slab  waveguides.  In  the  initial  stages  of  the  program,  numerous 
slab  waveguides  were  fabricated  in  LiNb03  y-cut,  z-propagating  substrates  by  standard 
Ti  (titanium)  indiffusion  techniques  [12].  Waveguide  quality  was  examined  by  incoupling 
experiments  utilizing  rutile  prism  coupling,  fiber  butt-coupling,  and  end-fire  coupling 
techniques.  In  all  cases,  the  coupling  was  successful,  and  low  scatter,  low  loss  waveguides 
were  achieved. 

The  second  critical  step  in  the  fabrication  of  a  parallel  waveguide  array  is  the  devel¬ 
opment  of  a  suitable  high  resolution  microlithography  process.  We  successfully  utilized 
a  liftoff  process  to  fabricate  high  resolution  channel  waveguides  in  LiNb03  from  a  set 
of  masks  obtained  from  Allied- Amphenol  Fiber  Optics  Products  Division.  One  mask 
contained  separate  3  pm,  5  pm,  7  pm,  and  9  pm  straight  waveguides  1  cm  in  length  in 
conjugation  with  a  two-port  directional  coupler  with  5  pm  channels  and  a  2  pm  gap. 
Each  of  the  channels,  as  well  as  the  directional  coupler,  were  shown  to  be  of  high  quality 
by  means  of  fiber  butt-coupling.  A  second  mask  contained  a  1  cm  long  Mach-Zender 
interferometer  with  5  pm  channels  and  a  3  pm  area  of  common  intersection,  which  was 
also  demonstrated  to  be  functional. 

Analysis  of  the  tradeoffs  between  rib  and  channel  waveguides  draws  one  to  the  con¬ 
clusion  that  although  channel  waveguide  arrays  are  potentially  much  easier  to  fabricate, 
rib  waveguides  may  prove  far  more  advantageous  for  our  application  due  to  the  low 
values  of  An  obtainable  from  Ti-indiffusion  in  LiNb03  and  the  correspondingly  large 
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between-channel  gaps  required  to  assure  effective  isolation.  Consequently,  we  ordered  a 
Commonwealth  Scientific  3  cm  ion  source  and  power  supply,  to  be  incorporated  into  our 
existing  Sloan  magnetron  sputtering  system  for  evaluation  of  the  feasibility  of  fabricating 
large  area  arrays  of  high  resolution  (8  //m  ribs  with  2  ^m  gaps)  rib  waveguides.  Such  an 
array  would  contain  1000  such  ribs  (range  channels)  in  1  cm  of  waveguide  cross-section. 

Finally,  in  support  of  our  efforts  to  upgrade  our  capacity  for  ion  beam  milling  ex¬ 
periments  with  in  situ  ion  beam  deposition,  ion-assisted  deposition,  and  planar  triode 
sputtering,  we  have  designed  (in  conjunction  with  L.M.  Simard  Co.,  Santa  Barbara,  C'A) 
and  purchased  a  multipurpose  vacuum  system  with  delivery  anticipated  in  19S7.  This 
unique  ultra-high-vacuum  system  will  allow  for  precision  ion  beam  milling  experiments 
with  a  fully  articulating,  cooled  substrate  fixture,  as  well  as  for  a  highly  flexible  array  of 
thin  film  deposition  and  processing  functions. 

D.  Waveguide  Outcoupling  Techniques 

Another  critical  issue  that  must  be  addressed  that  is  related  to  the  area  of  wave¬ 
guide  design  and  fabrication  is  the  development  of  a  technique  for  uniform  extraction  of 
light  from  each  waveguide  in  a  direction  normal  to  the  surface  of  the  integrated  optics 
substrate.  A  number  of  potential  methods  are  worthy  of  further  investigation,  including: 

1.  Uniform  Overlayer 

This  is  the  most  straightforward  method  for  generating  uniform  leakage  along  a  given 
waveguide.  In  this  case,  a  thin  film  overlayer  is  chosen  to  have  a  specific  refractive  index 
that  results  in  a  very  low  loss  per  unit  length  from  evanescent  wave  coupling.  If  the 
total  waveguide  loss  is  set  to  5%,  for  example,  then  an  essentially  uniform  vertical  light 
leakage  will  result.  This  mode  of  operation  is  inefficient  in  its  utilization  of  light  from 
the  laser  diode  source,  however,  since  95%  of  the  light  will  be  emitted  from  the  far  end 
of  each  waveguide  and  hence  will  not  contribute  to  the  detected  intensity. 

2.  Nonuniform  Overlayer 

If  the  refractive  index  of  the  waveguide  overlayer  is  varied  along  the  length  of  the 
waveguide,  the  loss  per  unit  length  can  be  increased  to  make  more  efficient  utilization  of 
the  incident  illumination.  The  advantage  of  this  method  is  the  gain  in  overall  through¬ 
put  efficiency,  but  potential  disadvantages  are  the  complication  of  achieving  a  specific 
refractive  index  profile  (or  equivalently,  a  thickness  profile  of  the  overlayer)  that  allows 
uniform  light  loss  per  unit  length,  and  the  difficulty  of  generating  sufficient  refractive 
index  variation  or  uniformity  of  thickness  control  to  achieve  the  desired  effect.  A  further 
inefficiency  that  complicates  the  uniform  overlayer  and  nonuniform  overlaver  methods 
is  the  isotropic  emission  into  a  full  hemisphere  of  the  evanescentlv  coupled  light  loss 


from  the  layer,  which  places  strict  demands  on  the  proximity  coupling  of  the  mask  and 
photodetector  required  to  capture  a  significant  fraction  of  the  emitted  light. 

3.  Uniform  Grating 

A  potential  avenue  of  inquiry  that  circumvents  a  number  of  the  limitations  of  the 
above  two  methods  involves  the  corrugation  of  each  waveguide  along  its  length  to  produce 
a  longitudinal  grating  coupler,  as  shown  in  Fig.  3.  In  this  case,  light  will  be  emitted 
from  the  waveguide  at  the  Bragg  angle  due  to  periodic  coherent  superposition  effects. 
This  allows  for  significantly  enhanced  latitude  in  detector  proximity  coupling,  and  can 
be  fabricated  with  high  uniformity.  If  the  modulation  depth  of  the  corrugation  is  fixed 
as  shown  in  Fig.  3,  the  modulation  depth  must  be  kept  small  to  maintain  a  uniform 
incident  intensity  down  the  length  of  the  waveguide. 

In  our  initial  experiments  with  grating  outcoupling,  we  have  utilized  direct  spin¬ 
coating  of  photoresist  on  top  of  slab  waveguides  in  LiNbC>3  to  produce  an  optically 
sensitive  outcoupling  layer.  Gratings  were  holographically  written  in  the  photoresist 
layer  by  means  of  a  specially  designed  interferometer  illuminated  by  a  He-Cd  laser  (441.6 
nm).  Grating  periods  of  1.4  pm  were  produced  routinely  in  this  arrangement. 

Successful  grating  outcoupling  was  obtained  over  large  areas  (4x6  mm)  with  this 
technique.  On  the  other  hand,  processing  parameters  required  for  obtaining  high  quality 
outcoupling  proved  extremely  difficult  to  stabilize.  This  is  principally  due  to  the  tradeoff 
between  the  modulation  depth  of  the  grating  (w’hich  must  be  low  to  afford  low'  efficiency 
but  uniform  coupling  down  the  entire  length  of  the  waveguide)  and  the  thickness  of  the 
photoresist  in  the  grating  valleys  (which  must  be  non-zero  to  promote  non-negligible 
evanescent  wave  outcoupling,  but  not  so  thick  as  to  exponentially  attenuate  the  evanes¬ 
cent  wave  amplitude  at  the  air-photoresist  interface).  In  the  follow-on  program,  we 
will  explore  ion  milling  of  the  grating  directly  into  the  waveguide  surface  as  a  means  of 
producing  a  more  reliable  grating  outcoupling  process. 

4.  Nonuniform  Grating 

Perhaps  the  most  promising  technique  of  all  involves  alteration  of  the  grating  modu¬ 
lation  depth  along  the  length  of  the  waveguide  to  allow  a  greater  fraction  of  the  light  to 
be  coupled  out  per  unit  length  (compensating  for  the  decreased  intensity  at  each  point  in 
the  waveguide  with  correspondingly  increased  modulation  depth).  Such  a  configuration 
is  shown  schematically  in  Fig.  4.  This  type  of  modulated  corrugation  could  be  produced 
by  holographic  interferometry  within  a  photoresist  layer,  with  longitudinal  variation  pro¬ 
vided  by  space-variant  reference  and  signal  beam  ratios.  The  photoresist  pattern  that 
results  would  then  be  subjected  to  ion  beam  milling  to  modulate  the  thickness  of  the  rib 
waveguide. 
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E.  Incorporation  of  the  Reference  Beam 

As  described  in  Appendices  A  and  B,  a  reference  beam  must  be  incorporated  within 
the  SAR  processor  architecture  to  re-encode  the  azimuth  phase  history  (doppler-shift)  as 
an  amplitude  modulation  suitable  for  detection  by  the  CCD  array.  During  the  contract 
period,  we  explored  several  potential  methods  for  incorporation  of  the  reference  beam. 
The  most  straightforward  of  these  is  simply  to  add  an  appropriate  frequency  RF  carrier 
to  the  SAR  return  signal,  which  will  produce  a  reference  beam  appropriately  doppler 
shifted  from  the  laser  frequency  so  that  it  can  interfere  with  each  separate  range-focused 
beam  at  the  entrance  of  the  parallel  rib  waveguide  array.  This  interference  will  produce 
the  requisite  amplitude  modulation  along  the  entire  length  of  the  array  provided  only 
that  the  coherence  length  of  the  pulsed  laser  is  of  order  the  length  of  the  array. 

F.  CCD  Detector  Interface 

Charge- coupled  device  arrays  with  adequate  performance  characteristics  are  com¬ 
mercially  available  at  the  present  time.  The  electronic  subsystem  necessary  to  convert  a 
commercial  CCD  array  (designed  as  a  closed-circuit  TV  camera)  so  that  it  can  operate 
in  the  shift-and-add  mode  required  by  the  SAR  processor  configuration  described  herein 
is  similar  in  nature  to  an  existing  design  fabricated  previously  at  Caltech.  The  key  issue 
here  is  the  interface  of  the  CCD  to  the  integrated  optical  chip.  Normally,  a  detector  array 
is  butt-coupled  at  one  end  of  an  integrated  optical  chip  (and  is  usually  a  one-dimensional 
detector).  In  the  proposed  processor  configuration  we  must  couple  the  light  out  from 
the  top  of  the  integrated  optical  waveguides  and  onto  a  two-dimensional  detector.  There 
are  (at  least)  two  possible  techniques  for  accomplishing  this.  The  first  is  by  means  of 
proximity  coupling.  In  this  method,  the  mask  and  CCD  combination  are  placed  as  close 
as  possible  to  the  surface  of  the  integrated  optical  rib  waveguide  structure.  The  CCD 
cannot  be  placed  in  actual  contact  with  the  10  chip  because  the  electrode  structure  of 
the  CCD  array  would  be  potentially  damaged,  and  in  addition,  the  characteristics  of 
the  optical  waveguides  might  be  affected.  One  possible  method  by  which  this  type  of 
coupling  could  be  effected  is  the  application  of  a  very  thin  (less  than  50  micron)  adhe¬ 
sive  layer  between  the  CCD  array  and  the  integrated  optical  chip.  This  layer  must  be 
nonconductive,  optically  transparent,  and  of  index  of  refraction  lower  than  that  of  either 
the  10  substrate  or  the  thin  film  overlayer  (if  such  a  layer  is  employed).  Alternatively, 
the  requisite  light  coupling  can  be  accomplished  by  means  of  a  thin  air  gap  produced  by 
microspacers. 

G.  Comparison  with  Present  and  Future  Electronic  Technology 

In  order  to  carefully  assess  the  potential  computational  throughput  gain,  as  well  as  the 
potential  weight,  power,  and  cost  advantages  of  the  IOSAR  processor,  we  also  initiated 
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during  the  contract  period  a  direct  comparison  between  current  (and  projected)  elec¬ 
tronic  processors  capable  of  synthetic  aperture  radar  image  formation  and  the  projected 
performance  of  a  fully  functional  (but  not  technologically  advanced)  IOSAR  processor 
architecture.  This  study  began  by  formulating  the  SAR  image  formation  problem  opti¬ 
mally  for  the  IBM  PC-XT  with  S0S7  math  co-processor  chip  as  well  as  for  the  YAX-VMS 
11/750  host.  Preliminary  results  indicate  a  tremendous  gain  in  overall  computational 
throughput  per  unit  power  for  the  IOSAR  architecture.  These  results  will  be  further 
documented  in  the  Final  Report  for  the  follow-on  contract  period.  The  comparative 
study  has  also  included  consideration  of  the  possibility  of  using  architecture  mapping 
algorithms  to  map  the  SAR  image  formation  problem  onto  a  single  YSLI  or  WSI  chip. 
These  efforts  have  not  met  with  success,  due  primarily  to  the  heavy  interconnection 
burden  required  by  the  image  formation  algorithms. 


H.  Extensions  of  the  Concept 


The  Integrated  Optical  Synthetic  Aperture  Radar  Processor  described  in  this  pro¬ 
posal  was  conceived  specifically  for  the  SAR  application.  However,  we  believe  that  the 
novel  concepts  that  we  have  discussed  herein  will  have  a  much  broader  significance. 
The  proposed  configuration  provides  for  the  processing  of  inherently  two-dimensional 
information  in  a  planar  integrated  optical  architecture.  The  information  is  processed  in 
both  the  transverse  and  longitudinal  dimensions  of  the  integrated  optical  chip,  and  de¬ 
tected  through  utilization  of  the  third  dimension  normal  to  the  plane.  This  innovation, 
coupled  with  the  time-and-space  integration  algorithm,  allows  extremely  complicated 
two-dimensional  operations  to  be  performed  at  very  high  potential  frame  rates.  There 
are  several  possible  variations  and  extensions  of  this  concept,  as  well  as  other  problems 
to  which  it  is  clearly  applicable. 


First,  methods  for  implementing  a  programmable  IOSAR  processor  are  of  interest. 
In  the  particular  implementation  drawn  schematically  in  Fig.  2,  the  azimuth  reference 
is  recorded  on  a  fixed  mask,  and  thus  cannot  be  reprogrammed  conveniently  for  SAR 
applications  involving  mobile  platform  geometries  and  dynamic  reference  conditions.  One 
such  method  is  the  replacement  of  the  fixed  mask  with  a  two-dimensional  spatial  light 
modulator  that  is  electronically  programmable  to  reflect  the  dynamic  flight  conditions 
requiring  compensation. 


Second,  the  direct  implementation  of  other  two-dimensional  optical  processing  oper¬ 
ations  can  be  achieved  utilizing  the  IOSARP  configuration  (with  modifications  in  some 
cases).  It  is  possible  to  implement  any  separable,  two-dimensional  signal  processing 
operation  using  this  approach.  For  example,  a  Fourier  transform  10  processor  for  image- 
based  inputs  is  conceivable.  The  implementation  of  non-separable  operations  such  as 
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two-dimensional  correlations  are  more  difficult,  but  still  quite  feasible.  In  particular, 
non-separable  functions  may  be  implemented  by  means  of  multiple  input  channel  10 
processors  modeled  after  the  single- channel  device  described  in  this  report. 
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Abstract.  A  real-time  optical  synthetic  aperture  radar  (SAR)  processor  is  de¬ 
scribed.  The  processor  utilizes  an  acousto-optic  device  as  the  input  electromc- 
to-optical  transducer  and  a  CCD  camera  that  serves  as  the  optical  detector  and 
simultaneously  performs  the  focusing  of  the  SAR  image  m  the  azimuth 
direction. 

Keywords.-  two-dimensional  signal  processing ;  synthetic  aperture  radar  acousto-optic 
devices;  imaging,  image  processing. 
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I.  INTRODUCTION 

Synthetic  aperture  radar  (SAR)  still  represents  the  most  successful 
application  of  optical  computing,  even  though  more  than  30  years 
have  elapsed  since  the  initial  demonstration  of  the  optical  SAR  pro¬ 
cessor.  Optical  computers  are  used  routinely  today  to  form  SAR 
images  from  radar  signals  that  are  collected  by  aircraft  or  spacecraft 
and  recorded  on  photographic  film.  The  requirements  of  modem 
radars,  however,  are  often  not  met  by  the  present  film-based  optical 
processors.  These  requirements  include  higher  image  quality  (resolu¬ 
tion.  dynamic  range,  artifacts,  etc.),  real-time  image  formation,  flexi¬ 
bility,  and  on-board  processing  capability.  In  order  to  meet  such 
requirements,  digital  computers  are  increasingly  being  used  to  form 
SAR  images.  The  dramatic  advances  in  microelectronics  in  recent 
years  have  made  it  feasible  to  construct  digital  SAR  processors  that 
can  provide  better  image  quality  (due  to  the  highest  accuracy  and  the 
flexibility  in  programming  a  digital  computer)  than  film-based  SAR 
processors.  In  addition,  digital  processors  can  be  built  (as  they  have  in 
some  cases)  that  have  real-time  and  on-board  processing  capabilities. 
It  is  not  the  purpose  of  this  paper  to  present  a  comparison  of  optical 
versus  digital  techniques  in  SAR  (the  interested  reader  is  referred  to 
Ausherman's  paper1),  but  it  is  necessary  to  point  out  that  a  number  of 
problems  remain  with  digital  SAR  processors.  The  complexity  of 
these  systems  results  in  a  very  costly,  power  consuming,  relatively  large 
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and  heavy  processor.  Thus,  even  though  it  is  conceivable  to  build  a  low 
resolution  real-time  digital  processor,  for  many  of  the  applications  this 
is  an  impractical  solution. 

Very  significant  advances  have  also  been  achieved  in  the  area  of 
electro-optics  which  have  resulted  in  corresponding  improvements 
in  the  state  of  the  art  of  optical  computers.  Specifically,  acousto¬ 
optic  technology  has  matured,  and  consequently  high  quality  de¬ 
vices  are  now  available  to  be  used  as  the  input  spatial  light  mod  ulator, 
semiconductor  detector  arrays  (CCD's  and  photodiode  arrays)  have 
proven  to  be  excellent  as  detectors  in  the  optical  computer,  and  the 
semiconductor  light  sources  (LED’s  and  laser  diodes)  which  have 
been  developed  primarily  for  optical  communications,  are  a  very 
efficient  light  source  in  many  optical  processors.  These  and  other 
developments  have  made  it  possible  to  construct  real-time,  high 
performance,  power  efficient,  compact  and  relatively  inexpensive 
optical  computers  such  as  the  acousto-optic  spectrum  analyzer1  and 
the  ambiguity  function  radar  processors.1  In  this  paper  we  describe 
an  optical  SAR  processor  in  which  the  modern  components  avail¬ 
able  for  optical  computing  are  fully  utilized,  and,  therefore,  we 
believe  that  the  resulting  processor  has  great  potential  for  applica¬ 
tions  where  on-board,  real-time  SAR  image  formation  is  required 

In  Sec.  2  the  fundamentals  of  SAR  are  reviewed,  in  order  to 
establish  the  notation  we  will  be  using  in  the  paper.  A  comprehensive 
treatment  of  SAR  can  be  found  in  Refs.  4  and  5.  In  Sec.  3  the 
operation  of  the  optical  processor  is  described  in  detail,  and  in  Sec.  4 
we  discuss  its  performance  characteristics.  Our  experimental  results 
to  date  are  reported  in  Sec.  5. 

2.  SAR  SIGNALS 

The  side-looking  SAR  geometry  is  depicted  in  Fig.  I.  The  vehicle 
(aircraft  or  spacecraft)  is  traveling  with  a  constant  velocity  v  along 
the  ^direction  and  at  height  h.  Coded  pulses  art  emitted  periodically 
through  an  antenna  mounted  on  the  side  of  the  vehicle.  In  this  paper 
we  assume  that  the  pulses  are  chirp  coded.  The  transmitted  wave¬ 
forms  S(t)  in  complex  notation  is  given  by 

S(t)  =  £  reel  exp[jb(t-nT)J]exp(j2irv0t)  .  (I) 

n 

where  r  is  the  pulse  duration.  I/T  is  the  pulse  repetition  frequency 
(PRF),  b  is  the  chirp  rate,  and  vQ  is  the  microwave  frequency  of  the 
radar.  The  transmitted  waveform  illuminates  a  patch  of  the  surface 
to  be  imaged  (the  ground),  and  pan  of  n  is  reflected  back  towards  the 
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vehicle  and  detected  by  the  antenna.  The  objective  of  S  AR  is  to  form 
an  image  of  the  ground  from  these  reflections.  We  will  discuss  how 
an  image  of  a  single  point  scatterer  (impulse)  located  at  coordinates 
f^and  ij  =  i)aon  the  ground  can  be  formed.  Since  the  S  A  R  image 
formation  is  a  linear  operation,  knowledge  of  the  impulse  response 
completely  describes  the  system.  The  signal  received  by  the  antenna 
due  to  a  single  point  scatterer  is 

r(t)  =  A(t)S(t-2R/c)  .  (2) 

where  R  is  the  instantaneous  range  from  the  antenna  to  the  point 
scatterer  at  (f0 .  »j_).  A(t)  is  the  far-fteld  pattern  of  the  antenna  From 
the  geometry  of  Fig.  I,  R  can  be  written  as 

R<«)  =  ✓<•»,- ■!)* +  R.1  =  +  RoJ  •  (3) 

where  R0  —  >/h^  +  (a2  is  the  range  when  the  vehicle  crosses  jj  =  rjQ. 
In  most  cases,  the  variation  of  the  range  R(t)  with  time  can  be 
neglected  within  the  duration  rof  a  single  radar  pulse,  and  therefore 
R(t)  3  R(nT).  In  addition,  the  range  R0  is  typically  much  greater 
than  the  distance 1 17  —  >70|.  Thu  allows  us  to  expand  the  square  root 
in  Eq.  (3)  as  follows: 

(vnT  -  .fo)2 

R(t)  =*  R(nT)  =  R0  + - -  (4) 

Using  the  above  equation  and  Eq.  ( I),  we  can  rewrite  Eq.  (2)  as 

[l  -  2R0/c  -  -  nT  I 

- - - J 

■  -1 

X  e*p[jb(t-2R0/c-ir-nT)J] 


Fip  2-  Acousto-optic/CCO  synthetic  epartur*  radar  procasaor.  (IS  — 
light  source.  AOD— acousto-optic  device.  F,  is  the  tocel  length  of  the  lens 
1^.  and  the  broken  lines  indicate  the  path  of  the  rays  of  the  reference 
wave  | 

r  t  -  2R0/c  -  nT  1 
r(t)  =  £  A(nT)  rect  - - - J 

n 

X  exp[jb(t— 2RQ/c  -  nT)2]  exp[-](2irv0/cR0)  (vnT-r^)2] 

X  exp(j2rrv0t)  (6) 

The  received  signal  r(t)  must  be  processed  to  produce  an  image  of  the 
point  scatterer  on  the  ground.  In  the  next  section  we  will  describe  in 
detail  how  the  real-time  optical  processor  that  is  presented  in  this 
paper  forms  the  SAR  image  from  the  signal  r(t).  In  deriving  Eq.  (6).  a 
number  of  simplifying  assumptions  and  approximations  were  made, 
which  may  or  may  not  be  justifiable  in  practice  depending  on  the 
particular  application.  We  will  proceed  with  the  description  of  the 
system  based  on  this  relatively  simple  form  of  the  SAR  equation  [Eq. 
(6)),  and  in  Sec.  4  we  will  discuss  the  necessary  modifications  to  the 
system  for  cases  in  which  the  conditions  set  in  this  section  are  not  met. 

3.  ACOUSTO-OPTIC/CCD  SAR  PROCESSOR 

A  schematic  diagram  of  the  processor  is  shown  in  Fig.  2.  The  system 
is  illuminated  by  a  pulsed  laser.  The  light  amplitude  of  the  source  is 
modulated  by 


X  exp(j2irv#t)exp(-j2irv0  2R0/c)exp(-j2rv0^r)  .  (5) 

The  variation  of  the  antenna  pattern  A(t)duringa  single  pulse  can  be 
neglected,  and  hence  A(t)  3  A(nT).  The  quantity  Sr  =  (vnT  - 
•?«?/  R0c  has  units  of  time  and  is  very  small  compared  to  2R0/c  if 
R«  »l  1  -  »»0I  Therefore,  Sr  can  be  neglected  in  Eq.($)  except  in 
the  term  exp(— j2rv0.ir)  because  in  (his  term  Sr  is  multiplied  by  v0 
which  is  the  very  Urge  radio  frequency.  The  term  exp<  -j4  rv*  RQ /  c) 
is  as  inconsequential  constant  phase  term  and  can  be  omitted.  With 
these  modifications,  we  can  write  the  final  form  for  r(t)  as 


P(t)  =  £  rect 
n 


(7) 


where  r0  is  the  duration  of  each  light  pulse.  The  spherical  lens  L, 
collimates  the  light  from  (he  source,  and  the  c>lindncal  lens  l.j 
focuses  the  light  in  the  vertical  (y)  direction  so  that  it  can  pass 
through  (he  aperture  of  the  acousto-optic  device  (AOD)  placed  in 
plane  P(.  The  received  radar  signal  r(t)  is  heterodyned  to  the  center 
frequency  v,  of  the  AOD,  and  a  reference  signal  Bexp(j2jrVjt)  is 
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added  to  it.  The  frequency  v2  is  chosen  such  that  the  difference  v2  — 
v,  is  equal  to  the  bandwidth  of  each  radar  pulse: 

v1  -  v,  =  b r/  v  .  (8) 

The  need  for  the  reference  signal  will  become  apparent  shortly.  The 
resulting  signal  r'(t)  is  related  to  r(t)  by  the  following  equation: 

r'(t)  =  r(t)  exp[-j2ir(va  -  v,)t]  +  Bexp(j2irv2t)  .  (9) 

r'(t)  is  applied  to  the  piezoelectric  transducer  of  the  AOD  of  Fig.  2. 
The  light  diffracted  by  the  AOD  is  modulated  by* 

S,(t.x)  =  rect  P(t)  r'  ^t  +  •  (10) 

where  x  is  along  the  direction  of  propagation  of  the  acoustic  wave,  X  is 
the  aperture  of  the  AOD,  and  va  is  the  speed  of  sound  in  the  AO  D.  The 
undiffracted  light  is  blocked  in  the  focal  plane  of  the  cylindrical  lens 
Lj.  The  combination  of  lenses  Lj  and  Lj  accomplishes  two  tasks.  First, 
since  the  two  lenses  are  separated  by  the  sum  of  (heir  focal  lengths,  a 
single  plane  wave  incident  on  Lj  will  be  recollimated  when  it  exits  the 
lens  Lj.  In  addition,  the  impulse  response  of  the  system  in  the  x 
direction  from  plane  P,  to  plane  P2  is  made  to  be  equal  to 

h(x,x)  =  exp£jb|(x ~ x)2]  exp(-j2rrv, x/va)  .  (II) 

where  x  is  the  horizontal  coordinate  in  plane  P2.  The  constant  bt  can 
be  set  by  appropriately  choosing  the  focal  lengths  of  L2  and  Ls  and 
the  distanced.  The  term  exp(— j2rrv,  x/v^)  reflects  the  fact  that  the 
optical  system  following  the  Bragg  cell  is  tilted  so  that  its  optical  axis 
coincides  with  a  wave  diffracted  at  the  Bragg  angle.  In  Eq.  (H)we 
assume  that  the  effects  of  the  lens  apertures  are  negligible.  In  the 
vertical  direction  y,  the  light  is  recollimated  by  the  lens  1.4,  and 
therefore  the  amplitude  of  light  entering  plane  P2  does  not  vary 
along  y.  At  P2  a  transparency  is  placed  immediately  in  front  of  the 
CCD  camera.  The  intensity  transmittance  T(x,y)  of  this  mask  is 


T(x.y)  =  -  +  -  cos  2»ru0x  + 


*2? 

— : —  • 

X 


where  u0  and  b^  are  constants  that  will  be  specified  later.  The 
instantaneous  intensity  incident  on  the  CCD  is  given  by 


I(x.y.t)  =  T(x. 


y)  |  f  si(‘- 


x)  h(x,x)  dx 


The  CCD  is  operated  in  the  shift-and-add  mode.7  In  this  mode  of 
operation  the  device  is  exposed  to  light  periodically.  The  photogen* 
erated  charge  is  accumulated  on  the  detector  for  a  short  time  inter¬ 
val.  and  after  each  exposure  the  entire  photogenerated  charge 
pattern  is  electronically  shifted  by  one  pixel  along  one  of  the  dimen¬ 
sions  of  the  CCD.  The  photogenerated  charge  due  to  the  next 
exposure  is  simply  added  to  the  charge  that  is  already  stored  in  each 
pixel,  and  then  the  process  is  repeated.  In  our  system  the  charge 
transferring  is  done  in  synchronism  with  the  PRF  of  the  radar,  and 
the  integration  penod  is  set  equal  to  the  period  T  of  the  radar.  The 
charge  generated  on  the  CCD  during  the  nth  radar  pulse  is  given  by 

t  +  aT/2 


‘.y.n)  =  J  i(*.y.t)  dt 


t-nT/I 


=  2Bv-rr#T(x,y)  A(nT)sinc 


27r(vj-v.)x  2rrv  ,  | 

X  cos  [  - -  +  (vnT  -  +  dj 

+  bias  terms  (14) 

The  derivation  of  Eq  ( 14)  is  straightforward  but  tedious,  and  it  is 
presented  in  Appendix  A.  This  charge  pattern  gets  shifted  along  the 
y  direction  by  (N  -  n)  pixels  on  the  CCD,  after  N  exposures  have 
taken  place.  The  total  charge  that  accumulates  at  each  pixel  located 
at  coordinates  (x,  y)  after  N  exposures  is 


Q(x.y)  =  £  Q  [x,y+(N-n)Ay.n] 
n 

Bv«rro  r br  /.  2R0va\ 

=  - r -  sine  - lx - I 

2  irva  \  C  / 

x  (y -  no  +  NAy)  cos  [< 


2Nv(JvJT? 


cos  [4>ru0x  +  <a] 


+  bias  terms  .  (15) 

In  the  above  equation  Ay  is  the  pixel  separation  in  the  CCD,  u0  =  (v2 
—  v, )/  va,  and  d  is  a  constant  phase  term.  Q(x,y)  is  the  final  output 
of  the  processor,  and  its  form  demonstrates  the  imaging  capability  of 
the  system.  In  the  x  direction,  Q(x.y)  is  a  sine  function  with  width 
rrv0/  brand  centered  at  x  =  2R0v0/c.  In  other  words,  we  obtain  the 
image  of  the  point  scatterer  located  at  f  =  fQon  the  ground,  focused 
in  the  f  direction.  The  resolution  in  the  (  direction  is  equal  to  St  = 
rrc/2br.  This  is  a  well-known  result  for  range  resolution  in  radar 
systems .*  The  carrier  at  spatial  frequency  2u0  arises  from  the  inclu¬ 
sion  of  the  reference  signal  in  the  acousto-optic  device  which  allowed 
us  to  record  on  the  CCD  the  phase  of  the  detected  signals  This  is 
necessary  in  this  system  because  the  Doppler  information  that  is 
essential  for  the  focusing  in  the  azimuth  (n)  direction  is  encoded  in 
the  phase  of  the  range  compressed  signal  Q(x.y,n).  In  addition, 
since  the  output  forms  on  a  earner,  it  can  be  easily  separated  from 
the  bias  term  in  Eq.  (IS). 

I  n  the  y  direction,  Q(x .  y)  is  also  a  sine  function  whose  position  is 
proportional  to  ij0.  the  location  of  the  point  scatterer  on  the  ground 
in  the  azimuth  direction  The  width  of  this  sine  function  is  equal  to 
cAyRg/NVjjt^T2,  and  it  determines  the  resolution  in  azimuth.  In 
ground  coordinates  rj,  the  width  of  the  sine  function  corresponds  to 
an  azimuth  resolution 
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D,  =  NvT  is  the  distance  that  the  vehicle  collecting  the  SAR  data 
travels  during  the  time  interval  NT,  and  it  is  equal  to  the  synthetic 
aperture  of  the  system.  Therefore,  the  azimuth  resolution  obtained 
with  this  processor  conforms  with  previously  derived  results’  for 
SAR  systems.  The  entire  pattern  in  shifted  in  the  y  dimension  by  the 
distance  NAy.  For  a  CCD  with  N  pixels  in  the  horizontal  direction, 
this  implies  that  after  the  signal  is  integrated  on  the  CCD  for  V 
pulses  to  produce  a  focused  image,  it  arrives  at  the  edge  of  the  device 
(at  y  =  NAy)  where  a  separate  CCD  stage  transfers  an  entire  line  of 
the  data  (a  slice  of  the  image  for  each  azimuth  position)  to  the  output 
pin  of  the  CCD  device  The  azimuth  slices  are  continuously  pro¬ 
duced  as  long  as  the  flight  continues,  producing  in  real  time  an  image 
of  a  long  strip  on  the  ground  parallel  to  the  direction  of  the  flight. 

4.  DISCUSSION 

In  this  section  we  discuss  some  of  the  issues  relative  to  the  perfor¬ 
mance  of  the  optical  processor  we  have  described.  We  will  address 
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several  performance  criteria  separately  and  comment  on  the  effects 
of  the  processor  architecture  and  the  devices  used  in  the  system  A 
complete  in-depth  analysis  of  the  performance  is  beyond  the  scope  of 
this  paper. 

4.1.  Aiimuth  resolution 

In  the  previous  section  we  derived  the  formula  for  the  resolution  that 
can  be  obtained  in  the  along  track  (azimuth)  direction: 


3.  =  7T - 7  '  (16) 

"*  2Nv0vT 

N  is  the  only  parameter  in  the  above  equation  that  is  related  to  the 
processor.  N  is  the  number  of  radar  pulses  that  are  used  to  form  the 
SAR  image  of  each  point  scatterer  on  the  ground.  In  the  AOf  CCD 
processor,  N  is  also  equal  to  the  number  of  pixels  of  the  CCD  in  the 
shifting  direction  (y).  Commercially  available  CCD  cameras  have 
approximately  500  pixels.  Several  devices  with  approximately  twice 
this  number  have  recently  been  developed*  or  are  currently  under 
development.  If  a  CCD  with  10*  pixels  is  used,  the  azimuth  resolu¬ 
tion  obtainable  with  this  processor  is  22.5  meters  on  the  ground.  This 
figure  was  obtained  by  using  the  parameters  of  the  SEASAT  radar4 
that  was  flown  by  NASA.  The  azimuth  resolution  will,  of  course,  be 
different  for  different  radars.  To  obtain  lower  resolution,  larger 
CCD’s  would  have  to  be  used.  It  is  unlikely  that  a  significant  increase 
in  the  size  of  monolithic  CCD  cameras  will  take  place  in  the  near 
future.  It  is  possible,  however,  to  construct  larger  arrays  by  inter¬ 
facing  a  number  of  separate  chips.  As  N  is  increased  and  finer 
resolution  is  achieved,  other  factors  in  the  processor  become  limit¬ 
ing,  such  as  optical  aberrations,  mechanical/ electronic  stability,  and 
light  source  coherence  over  the  lime  period  NT.  Finally,  for  low 
resolution  imaging,  some  of  the  approximations  that  were  made  in 
Sec.  2  may  not  be  valid  For  instance,  if  the  approximation  that  the 
range  of  each  point  target  remains  constant  as  the  plane  flies  past  it  is 
not  valid,  an  effect  known  as  range  curvature*  results.  Range  curva¬ 
ture  can  be  corrected  in  our  processor  simply  by  tilting  one  of  the 
cylindrical  lenses  in  the  system  which  introduces  barrel  distortion 
that  can  exactly  cancel  out  the  range  curvature. 

4J.  Image  size  In  azimuth 

The  AO/ CCD  processor  continuously  produces  slices  of  the  SAR 
images  in  the  azimuth  direction,  and  therefore  the  size  of  the  SAR 
image  is  not  limited  by  the  processor  in  this  direction. 

4J.  Range  resolution 

The  resolution  in  the  range  direction  (f)  on  the  ground  is  given  by 


=  ^  (,7> 

The  bandwidth  of  the  chirped  radar  pulses  is  equal  to  br/  tr,  and 
therefore  the  resolution  in  range  is  limited  by  the  temporal  band¬ 
width  of  the  processor.  The  bandwidth  of  the  AOD  used  in  the 
system  is  not  the  primarily  limiting  factor  since  devices  with  several 
gigahertz  bandwidth  are  available.  More  importantly,  the  duration 
of  the  light  pulse  r0  must  be  shorter  than  the  inverse  of  the  band¬ 
width  of  the  signal  to  avoid  smearing  of  the  range  compressed  signal 
(Eq.  ( 14)}.  For  typical  radar  bandwidths,  r0  should  be  several  nano¬ 
seconds  .  Several  laser  sources  are  capable  of  producing  nanosecond 
pulses  with  adequate  peak  power.  We  are  primarily  considering  laser 
diodes  because  they  are  compact  and  power  efficient.  Laser  diodes 
can  las*  in  a  single  mode  under  cw  operation  (and  therefore  be 
coherent),  but  when  they  are  pulsed,  several  modes  exhibit  gain 
during  the  transients.  The  need  for  coherence  in  the  system  imposes 
with  present  sute-of-the-en  devices  a  minimum  pulse  width  of  —20 
m.  Therefore,  if  a  laser  diode  is  used  in  the  system,  the  bandwidth  of 
the  radar  must  be  equal  fo  50  MHz.  which  corresponds  to  range 
tesolution  equal  to  3  m.  This  can  be  lowered  by  an  order  of  magni¬ 


tude  if  a  pulsed  gas  !a>er  ri  u->ed  The  additional  factors  mentioned  in 
our  discussion  on  azimuth  resolution  also  become  increasingly 
important  for  finer  range  resolution. 

4.4.  Image  size  in  range 

The  number  of  pi  sets  in  the  range  direction  is  limited  by  the  space- 
bandwidth  product  of  the  AOD  or  the  CCD  in  the  nonshifting 
direction,  whichever  is  smaller.  At  the  present  time,  the  CCD 
imposes  this  limit  at  approximately  ICP  pixels.  The  range  coverage 
can  be  extended  by  using  several  processors  operating  in  parallel. 

4.5.  Flexibility 

A  potential  application  for  the  AO;  CCD  processor  is  on-board 
real-lime  image  formation  In  such  applications  the  radar  geometry 
can  be  expected  to  change  from  time  to  time  or  even  continuously. 
Enough  flexibility  must  be  built  into  the  system  so  that  it  can  be 
adaptable  to  the  changing  conditions.  One  method  for  achiev  mg  this 
is  by  the  use  of  a  real-time  2-D  spatial  light  modulator*  in  place  of  the 
fixed  mask  in  Fig.  2.  The  impulse  responseof  the  processorcan  then 
be  altered  to  track  changes  in  the  radar  geometry. 

4.6.  Dynamic  range 

The  dynamic  range  is  defined  as  the  maximum  unsaturated  bright¬ 
ness  in  the  SAR  image  divided  by  the  minimum  detectable  level 
above  the  noise.  In  general,  a  large  dynamic  range  can  be  observed  in 
a  SAR  image  because  the  noise  detected  at  radio  frequencies  is  very 
effectively  suppressed  through  the  pulse  compression  in  both  range 
and  azimuth  that  is  performed  to  bring  the  SAR  image  into  focus. 
The  potential  for  high  dynamic  range,  however,  can  only  be  realized 
if  the  processor  that  forms  the  SAR  image  can  provide  it.  In  an 
optical  processor,  essentially  every  component  of  the  sy  stem  con¬ 
tributes  to  diminish  the  dynamic  range  available  at  the  output  by- 
adding  noise  or  by  introducing  nonlinearities  and  aberrations  in  the 
system  In  the  AO/ CCD  processor,  the  CCD  detector  is  the  domi¬ 
nant  factor  that  limits  the  dynamic  range.  The  maximum  bright  ness 
of  the  SAR  image  is  equal  icvthe  maximum  charge  that  can  be  stored 
in  each  pixel  of  the  CCD,  minus  the  charge  that  corresponds  to  the 
bias  terms  in  Eq.  (15)  and  the  dark  current  that  accumulates  during 
the  integration  time.  The  bias  is  signal  dependent  and  is  largest  at 
portions  of  the  SAR  image  that  have  maximum  average  intensity. 
Therefore,  relatively  sparse  objects  can  be  imaged  with  a  larger 
dynamic  range.  For  a  single  point  scatterer  on  the  ground,  the  bias 
can  be  only  1/2  of  the  full  dynamic  range  of  the  CCD,  but  for  a 
distributed  target  the  bias  increases  as  v/M,  where  M  is  the  average 
number  of  point  scatterers  on  the  ground  that  are  located  at  the  same 
range  in  a  distance  NvT  along  the  azimuth  direction.  The  minimum 
detectable  level  at  the  output  of  the  processor  is  determined  by  the 
fixed  pattern  noise  of  the  CCD  (the  variation  of  the  dark  current 
from  pixel  to  pixel),  optical  scattering,  quantum  noise  primarily  due 
to  the  bias,  and  thermal  and  amplifier  noise  from  the  electronics.10 
Our  calculations"  show  that  the  dynamic  range  of  the  images  that 
will  be  produced  by  the  experimental  processor  we  are  presently 
constructing  is  —200  I.  A  commercially  available  CCD  camera  is 
used  in  the  experimental  processor,  and  post-detection  signal  pro¬ 
cessing  is  applied  in  order  to  minimize  noise.  The  dynamic  range  can 
be  increased  to  over  10*:  I  by  cooling  the  CCD  to  eliminate  the  dark 
current  and  by  using  a  CCD  specifically  designed  for  this  application. 

5.  EXPERIMENT 

The  CCD  camera  is  utilized  in  this  system  not  only  as  the  optical 
detector  but  also  as  a  multichannel  correlator  by  operating  it  in  the 
shift-and-add  mode  described  earlier.  The  signal  that  comes  out  at 
each  column  (at  position  x)  of  the  device  is  the  correlation  between 
the  transmittance  of  the  mask  as  a  function  of  y  at  the  same  position 
x  and  the  signal  corresponding  to  the  temporal  variation  of  the 
incident  light  intensity.  We  have  performed  an  experiment  to  dem¬ 
onstrate  the  operation  of  the  CCD  as  an  optical  multichannel  corre¬ 
lator.  The  experimental  setup  is  shown  in  Fig.  3  A  light  emitting 
diode  (LED)  is  used  as  the  light  source.  The  intensity  of  the  LED  is 


Fig.  3.  Eiparimantal  alac'ro  optic  ambiguity  (unction  procasio r. 
(LED— light amitting  diode.  L,  —collimating  lam.  M — mash,  Lj  —imaging 
Ians.  CCD— camara.) 
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Fig.  4.  Ambiguity  processor  signals:  (a)  input  (LEO)  wavaform;  (b)  mask; 
(c)  ambiguity  function. 

temporally  modulated  by  the  chirp  signal  shown  in  Fig.  4(a).  The 
chirp  signal  that  modulates  the  LEO  represents  the  variation  of  the 
intensity  detected  by  the  CCD  (the  function  Q(  x .  y .  n)  in  Eq.  ( I 4)]  as 
a  function  of  pulse  number  n.  The  duration  of  this  waveform  is  equal 
to  NT.  and  us  bandwidth  is  equal  to  v0vJNT/cR0  1  The  bandwidth 
of  the  waveform  in  Fig.  4(a)  is  30  Hz.  and  its  duration  is  I  s.  The 
CCD  we  used  in  this  experiment  has  N  =  512  pixels  in  the  vertical 
dimension,  and  therefore  the  duration  of  one  second  corresponds  to 
a  radar  pulse  repetition  frequency  equal  to  l/T  =  05  kHz.  These 
parameters  were  chosen  to  be  approximately  equal  to  the  parame¬ 
ters  of  an  experimental  aircraft  radar  which  is  flown  periodically  by 
the  Jet  Propulsion  Laboratory.  We  intend  to  use  data  collected  by 
this  radar  to  test  the  processor  when  it  is  completed. 

The  light  from  the  LED  is  collimated  by  lens  L,  in  Fig.  3,  and  it 
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illuminates  the  mask  denoted  by  M  in  Fig  J.  A  photograph  of  this 
mask  is  shown  in  Fig.  4(b).  The  intensity  transmittance  in  the  vertical 
direction  is  a  chirp  signal  whose  chirp  rate  is  proportional  to  the 
horizontal  position.  The  signal  at  the  center  of  the  mask  matches  the 
waveform  that  is  used  to  modulate  the  LED.  The  mask  is  imaged  by 
lens  Lj  onto  the  CCD.  which  is  operated  in  the  shifi-and-add  mode 
The  output  signal  from  the  CCD  at  each  horizontal  position  is  the 
correlation  of  the  waveform  modulating  the  LED  and  the  transmit¬ 
tance  of  the  mask  at  that  horizontal  position.  The  output  waveforms 
from  the  CCD  from  selected  vertical  locations  are  shown  in  Fig  4(c). 
This  2-D  waveform  is  known  as  the  ambiguity  function11  of  the  chirp 
signal,  and  Fig  4(c)  displays  the  familiar  tilted-ridge  structure  of  the 
ambiguity  function  for  this  particular  waveform. 
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7.  APPENDIX  A:  DERIVATION  OF  EQUATION  (14) 


t  +  nT/2 

Q(x.y.n)  =  f  l(x.y.t)  dt  . 


(14) 


t-nT/2 

Substituting  Eq.  ( 1 3)  in  Eq.  (14): 

l  +  nT/2 


QU 


,y.n)  =  T(x,y)  J  ^S,(t.x)  h(x.x)  dx 


t-nt/2 


dt  .  (Al) 


We  will  first  evaluate  the  integral  over  x  in  the  above  expression. 
Using  Eqs.  (10)  and  (II): 


g(t,x)  =/s((t.x)h(x.i)dx 

-jrtet  )  P(t)P  (i  +  exp(jbt(x  -  x)J] 

/  2rrv,x\ 

X  exp  ^-j  — - - jdx  .  (A2) 


Substituting  the  expressions  for  P(t)  and  Pfi)  from  Eqs.  (7)  and  (9): 


si,.;)  =  x  ■“>  f'm  (?) 

n 

H  (,+^)]} 

X  exp[jbl(x-x)i]  exp  dx  .  (A3) 


+  Bexp 


and  the  expression  for  the  received  signal  r(t)  from  Eq  (6)  is 
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/  l  —  nT  \  r  /  x  \  (  ft  -  2R  c  -  nT  +•  x.  v 

=  Zrect  [ - ~ - )jrta  lx/l2;A<nT)  ret‘  [ - - - J«p[jb,t  -2R0  c  nT  fu,!1] 

[2rr\Q  1 

_J  cR  <vnT  -10n{*P  [j2»rv0(t  +•  *  vo)]  exp  [  ~)2rr(v0  -  v,)  (I  +  X  vo»]  4-  Bexp  (j2tv,(1  -I-  x/vB)] 
X  exp  [jb,<x  -x)2]  exp(-j2rrv,x/va)  dx 


The  laser  pulse  width  r0  is  chosen  sufficiently  smaller  than  the  bandwidth  of  the  radar  pulses  <  rQ  <  br  ir)  so  that  the  temporal  sanation  of  the 
signal  in  the  above  equation  can  be  neglected  over  the  duration  of  each  light  pulse  With  this  approximation  Eq  ( A4)can  be  rewritten  as  folio »> 


t.x)  =  J  rect  ^ - - - j  exp(j2 rrv, t)  j rect  (  — )|v  A(nT)  rect  £ - — —  j  exp  [jblx  va  -  2R0  c)2] 

*  exP  (vnT  -  rjj,)2^  4-  Bexp  [j2»r(v2  -  v()x,  vff]  exp  [;2  rr<  v2  -  v, )  nT]  j  exp  [jb,  ( x  -  x); ]  dx  . 


We  can  assume  that  |  X|  >  var  4-  2R0va/c,  i  e  .  the  range  delay  is  such  that  the  reflected  radar  pulse  is  entirely  in  the  aperture  of  the  AOD  at  the 
instant  the  laser  is  pulsed.  This  can  be  ensured  by  appropriately  setting  the  relative  phase  of  the  radar  and  laser  pulse  repetition  frequencies  1  he 
parameter  6,  (the  focusing  power  of  the  optical  system  in  the  x  direction)  is  set  equal  to  b/  vo2  Finally,  the  offset  frequency  v,  —  v,  is  made  to  be 
equal  to  nv  T,  where  m  is  an  integer  This  can  be  accomplished  by  deriving  v,  -  v2  from  the  pulse  repetition  frequency  I ;  T  The  function  g(x  .t) 
can  then  be  written  as  follows: 


g(x.  t)  -  V  rect 


(~ ~)jA,nT)  exp  exp(J“7^)/r 


*''vo-:R0'c 


X  exp  -J - J- (x -2R0va,  c)  x  dx  +  Bexp  [j2rr(v2  -  v,)x,  v„]  exp  [jn^fvj  -  v,)2,  b]>  = 


X  \  A(nT)  varsinc 


[— ,i-2R°Vc)]  exp^-j  exp  ~  {\ -  2R0v„/c)  (2R0va,  c)j 


X  exp  -j 


~J7T(*J-4RoJv<»J'ci)  +  Bexp[J2rr<v)  -v2)x/vj  exp  [j  r^fVj  -  v, b]  j 


We  substitute  the  above  expression  in  Eq.  ( 14): 
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[br  .  1  r2rr(vj-v,)i  2rrv0  l) 

—  (*-2R0v./c)|co»  I - — -  +  (vnT  -  i\Qf  4  ^  I  >  Jt  =  2T(i.»Bvorr0A(nT) 

*,nc  ,*_2R°V®/C)j  co*  [2®fy2  _v|)»  Ya  +  2rrv0(vnT  -  n,,)2  cR„  4-$]  4-  bias  term* 
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In  the  abo \e  equation,  the  phase  terms  from  Eq.  (A6)  that  do  not 
depend  on  the  index  n  have  been  lumped  in  the  parameter 

8.  APPENDIX  B:  DERIVATION  OF  EQUATION  (15) 


Q(*.>)  =  2  Q(x.y-nAy  +  NAy,nT)  .  (15) 

■ 

Substitute  Eq.  (14)  in  the  above  equation  and  the  expression  for  the 
transmittance  of  the  mask  from  Eq.  (12): 

*  (II 

Q(*.y)  =  X  2Bv-rr«  \i  +  I 

m 

X  cos  jrbj(y —(n  — N)  Ay)2/x  J  j 

X  A(nT)  sine  ■  (x-2R0va/c)J 


X  cos[2rr(vJ-v()Vv0  +  2Bnr0(vnT-i70)1/cR(,  +  d>3  .  (Bl) 

The  spatial  frequency  u0  on  the  mask  must  be  selected  so  that  the 
sum  frequency  u#+(vj  —  v,)/va  of  the  two  cosine  functions  in  the 
above  equation  is  sufficiently  large  so  that  the  signal  term  that  is 
riding  on  the  sum  frequency  can  be  separated  from  the  bias  terms  in 
the  spatial  frequency  domain.  To  accomplish  this,  we  set  u0  = 
(Vj  —  v,)/ va.  We  write  the  product  of  the  cosines  as  a  sum,  and  the 
term  of  interest  in  Eq.  (Bl)  takes  the  following  form: 

Q(i^y)  38  (Bvarr0/2)sinc[(br/jrv0)  (x-2R0va/c)] 

N 

X  X  A(n"0  co*  [4«xt0x—  irbjjf2/! 

a 

+  2«-bj(n—N)Ayy/x— irbjA^fn  — N'),/i  +  2irv0vJTJnJ/cR0 

-  4«r*0unTn0cR0  +  2rv0r,0,/cR0  +  *]  .  (B2) 

The  focusing  of  the  SAR  image  in  the  azimuth  direction  (ij)  is 
accomplished  by  the  summation  over  n  in  the  above  equation.  The 
summation  will  result  in  a  well-focused  image  in  the  y  direction  if  the 
terms  containing  n2  in  the  argument  of  the  cosine  are  canceled.  To 
accomplish  this,  the  constant  bj  of  the  mask  is  set  appropriately. 

(2»#vJT,/eRa)nI  =  (bjAyJ/xJn1  — 


-b,  = 


Zv^T1* 


2v0v2T> 

“75 — <2RoVc> 


4v#vJT,vi 


Equation  (B3)  and  the  value  we  previously  chose  for  u0  provide  all 
the  information  necessary  to  construct 'the  mask  that  is  used  in  the 
system.  Equation  (B3)  is  substituted  in  Eq.  (B2),  and  all  inconse¬ 
quential  phase  terms  are  included  in  the  parameter  $.  Equation  (B2) 
then  reduces  to  the  following  equation: 

Q(x.>)  =  (Buorr0/2)  sinc[br/»rva(x-2R0va/c)] 

N 

X  X  A(nT)  cos  [4jru0x  +(4irv0vJTJ/ciyR0) 


X  (y-(-70Ay/vT)  +  NAy)n+*}  (B4) 

The  summation  over  n  in  Eq.(B4)  is  recognized  as  the  discrete  cosine 
transform^oftheantenna pattern  A(nT).  In  ordertogetan  approx¬ 
imate  result,  we  can  assume  that  the  antenna  pattern  is  so  broad  in 
the  azimuth  direction  that  its  main  lobe  uniformly  illuminates  a 
distance  vNT  on  the  ground.  Therefore,  A(nT)  is  approximately  a 
constant  forO>n>N,  and  Eq.  (B4)  can  be  evaluated  as  follows11: 


Q(*.y)  =  (BvarrQ/2)  sine  [(br/ rrvj  (x-2R0vo/c)} 

sin  [UNm^T^/cAyR,,)  (y-(*r0Ay/vT)  +  NAy)] 
sin  [(2irv{}v2T1/cAyR0)  (y-(»j0Ay/vT)  +  N  Ay)} 

X  cos  (4iru  x  +  <$)  .  (E 


The  ratio  of  the  sines  in  Eq.  (B5)  can  be  approximated  by  a  sine 
function  for  large  N.  and  this  leads  us  to  the  final  form  of  Eq.  (15): 

Q(*.y)  =  (NBvarr0/2)sinc[(br/irva)  (x-2R<iv„/c)3 

X  sine  [(2NvQvJTJ/cAyR0)  (y-(n0Ay/vT)  +N Ay)} 

X  cos  (4mi0x  +  £)  .  (15) 
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ABSTRACT 


Real  time  acousto-optic  SAR  processors  are  described  and  experimentally 
demonstrated.  SAR  imaging  is  performed  in  one  of  the  architectures  by  applying 
the  signal  to  an  acousto-optic  device  and  corelating  it  with  chirp  signals  recorded 
on  an  optical  transparency  by  time  integration  on  a  CCD  detector.  In  a  different 
implementation,  the  imaging  is  performed  by  interfering  the  light  beams  diffracted 
from  two  separate  acousto-optic  devices,  one  modulated  by  the  radar  signal  and  the 
second  by  the  reference  chirp  waveform. 


INTRODUCTION 


Synthetic  Aperture  Radar  [l,2]  is  a  powerful  technique  for  the  formation  of 
high  resolution  images  at  long  wavelengths.  The  most  attractive  feature  of  the 
technique  is  the  fact  that  the  resolution  is  not  limited  by  the  physical  aperture 
of  the  receiving  antenna,  but  the  bandwidth  of  the  radar  and  the  distance  the  vehicle 
carrying  the  antenna  travels  while  each  point  of  the  object  being  imaged  remains 
illuminated  by  the  radar.  Therefore,  SAR  images  of  the  earth  can  be  formed  with 
excellent  resolution  by  aircraft  or  spacecraft  that  can  carry  only  small  antennas. 
Until  recently,  the  usefullness  of  the  technique  has  been  limited  to  applications 
that  do  not  require  real  time  image  formation  because  of  the  high  computational  load 
associated  with  forming  the  focused  image  from  the  received  radar  signals.  Typically, 
the  radar  returns  are  recorded  on  the  aircraft  or  spacecraft  carrying  the  radar  and 
the  images  are  formed  later  on  the  ground  by  optical  or  digital  computers.  The 
requirements  of  modern  radars  and  advances  in  signal  processing  technologies  have 
raised  the  possibility  of  producing  SAR  images  in  real  time  on-board  the  vehicle 
carrying  the  radar.  Such  capability  would  enhance  significantly  the  applicability  of 
the  SAR  imaging  method.  The  advances  in  the  state-of-the-art  of  digital  and  optical 
signal  processing  technologies  that  have  been  accomplished  in  recent  years  make  it 
possible  to  form  SAR  images  in  real  time.  Processing  speed  however  is  only  one  of 
the  requirements  for  a  practical  SAR  processor.  Factors  such  as  power  consumption, 
weight,  size,  cost  and  flexibility  are  also  very  Important.  The  optical  implemen¬ 
tation  is  at  a  relative  disadvantage  to  a  digital  Implementation  with  respect  to 
flexibility  and  accuracy,  but  low  power  consumption,  size/weight  and  cost  make  the 
optical  implementation  the  only  practical  solution  in  many  SAR  imaging  applications. 
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We  have  described  recently  [3,4]  an  optical  architecture  for  SAR  imaging  that 
is  implemented  with  an  acousto-optic  device  (AOD)  as  the  input  transducer  and  a  2-D 
CCD  detector  array  on  which  the  SAR  image  is  formed.  In  this  paper  the  principles  of 
operation  of  this  architecture  are  reviewed  and  experimental  demonstration  of  real 
time  SAR  imaging  is  presented.  A  new  interferometric  architecture  Is  also  prsented 
in  which  the  Impulse  response  of  the  system  is  controlled  with  a  second  AOD  that  is 
modulated  by  a  chirp  waveform.  The  new  architecture  provides  the  flexibility  of 
dynamically  adapting  the  parameters  of  the  optical  processing  to  track  changes  in 
the  radar  geometry. 


EXPERIMENTAL  PROCESSOR 


We  will  present  the  operation  of  the  processor  by  describing  how  it  produces  a 
focused  image  of  a  point  target.  A  focused  SAR  image  is  related  to  the  reflectivity 
of  the  area  being  Imaged  through  a  2-D  linear  operation.  Therefore  the  impulse 
response  describes  the  system  completely.  The  typical  side-looking  SAR  geometry  is 
shown  in  Figure  1.  An  airplane  flying  at  a  constant  velocity  vfl  parallel  to  the 

ground,  illuminates  an  area  on  the  ground  through  a  small  antenna  mounted  on  the  side 
of  the  aircraft.  The  illuminating  radiation  is  temporally  modulated  by  periodic 
linear  FM  pulses.  The  transmitted  waveform  is  given  by  the  following  expression. 

s(t')  »  I  rect  R  -  exp  C  jb1  (t'-nT)2]  exp  ( jcjQt  * ) .  (1) 

n  o 


t'  is  time,  1/T  is  the  pulse  repetition  frequency  (PRF),  tq  is  the  duration  of  each 
pulse,  b^  is  the  chirp  rate  and  wq/2tt  is  the  center  frequency  of  the  radar.  The 


single  point  scatterer  on  the  ground  scatters  the  illuminating  field  and  part  of  it 
is  detected  by  the  antenna  on  the  aircraft.  For  a  single  point  scatterer  the 
received  waveform  is  a  replica  of  what  was  transmitted  delayed  by  2R/c,  where  R  is 
the  distance  from  the  antenna  to  the  scatterer  and  c  is  the  speed  of  light.  R 
changes  with  time  as  the  aircraft  flies  by  the  target  and  from  Fig.  1  we  find  that 

(l+vaJ  (f  -  to)2/  2K02). 


R(t’)  - 
where  R 


2  2 
<R  +  x  ) 
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*  -  [R 


2  2 
+  v 

o  a 


is  the  distance  to  the  scatter  at  the  time  t' 


t  when  the  vehicle  crosses 
o 


over  the  scatterer  and  the  assumption  has  been  made  that  Rq  >>x. 


Typically  the 


variation  of  R(t')  during  the  time  the  scatterer  is  illuminated  is  small  enough  and 
it  does  not  affect  the  time  of  arrival  or  the  chirping  rate  of  the  received  pulses. 
Therefore  each  pulse  is  received  with  approximately  the  same  delay,  2R  /c.  This 


delay  is  measured  by  compressing  the  received  pulses  and  the  image  of  the  scatterer 
is  focused  in  this  manner  in  the  across  track  (or  range)  dimension.  The  variation 
of  R  appears  as  a  phase  shift  on  the  received  waveform.  This  phase  modulation  is 
then  used  to  focus  the  image  in  the  along  tract  (or  azimuth)  direction.  With  these 
assumptions  and  ignoring  the  effects  of  the  antenna  pattern,  the  received  waveform 
can  be  written  as  follows. 
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( t*  -  2R  /c  -  nT) 

s (t *  -  R(t'))  =  Z  rect  - 2 -  exp  (t*  _  2Rq/c  -  nT)] 

n  L  o  J 

x  exp  (Joj  t’)  exp  [jw  v  2  (nT  -  t  )2/R  ].  (2) 

u  os  o  o 

It  Is  customary  to  think  of  this  long  1-D  signal  as  a  2-D  unfocused  image  of  the 
ground.  We  define  x  ■  t‘-nT  and  t  *  nT  and  we  write  the  received  signal  as  a  2-D 
function  as  follows. 

[T  -  2R  /cl  - 

— - 2 —  exp  [jt^  (t  -  2Rq/c)  ] 

x  exp  [Joj  v  2  (t  -  t  )2/R  ]  exp  [joi  x].  (3) 

o  a  o  o  o 

The  above  2-D  function  is  in  general  a  non-symmetric  zone  plate  whose  center 
(t  ■  2Rq/c,  t  *  tQ)  is  proportional  to  the  position  of  the  point  scatterer  on  the 

ground.  In  the  optical  processor  to  be  described  this  signal  is  focused  by  spatial 
integration  in  the  range  (t)  direction.  The  variable  t  does  not  represent  the 
continuous  time.  It  is  a  discrete  time  variable  sampled  at  the  radar  PRF.  We  will 
refer  to  it  as  simply  time  in  the  rest  of  our  discussion.  In  order  to  characterize 
and  demonstrate  the  operation  of  the  real  time  SAR  processor,  we  have  built  an 
apparatus  that  generates  the  Impulse  response  signal  of  Eq.(3).  A  schematic  diagram  of 
the  impulse  response  generator  is  shown  in  Fig.  2.  The  frequency  of  a  10  MHz 
oscillator  is  divided  by  a  factor  of  2x10^  to  produce  a  phase  locked  PRF  clock  at 
500  Hz.  The  same  10  MHz  frequency  is  used  to  derive  a  60  MHz  signal  and  this  signal 
is  then  applied  to  a  digitally  controlled  phase  shifter.  The  phase  shifter  is 
addressed  by  a  read-only  memory  circuit  in  which  a  quadratic  function  is  stored.  The 
quadratic  function  stored  in  the  memory  is  displayed  in  Fig.  3.  This  was  obtained 
through  a  D/A  converter.  The  memory  chip  is  read-out  at  the  PRF  of  the  radar  and 
consequently  the  signal  following  the  phase  shifter  has  a  slow  linear  frequency 
modulation  (20  Hz/frame  in  the  experiments  presented  later)  simulating  the  t  varia¬ 
tion  of  the  signal  in  Eq.  (3).  The  PRF  clock  triggers  a  circuit  that  produces  a 
very  short  pulse  (16. -6  nseconds)  which  is  in  turn  mixed  with  the  phase  shifted 
60  MHz  carrier  and  applied  to  a  surface  acoustic  wave  device.  The  dispersive 
properties  of  the  SAW  device  cause  the  impulses  to  be  spread  in  time  to  6  useconds 
and  in  frequency  to  the  range  50  to  70  MHz.  The  SAW  introduces  in  this  manner  the  x 
modulation  of  the  signal  in  Eq.  3.  The  output  of  the  SAW  is  a  replica  of  the  signal 
of  the  signal  that  Is  received  by  a  synthetic  aperture  radar  illuminating  a  scene 
containing  only  one  point  reflector.  This  signal  is  displayed  by  mixing  it  to  base¬ 
band  and  modulating  the  intensity  (z  axis)  of  an  oscilloscope.  The  horizontal  (x) 
amplifier  of  the  oscilloscope  is  triggered  by  the  PRF  clock  and  the  vertical  axis  is 
deflected  by  a  slow  ramp.  Fig*  4  shows  the  2-D  zone  plate  produced  by  the  system 
described  above.  This  zone  plate  is  a  representation  of  the  unfocused  image  of  a 
point  scatter. 

The  signal  from  the  impulse  response  generator  is  applied  to  an  AOD.  A  shear 
wave  TeO^  device  was  used  in  the  experiments  with  70  ysec  delay  window  and  30  MHz 

(40-70  MHz)  3dB  bandwidth.  The  AOD  is  illuminated  by  a  pulsed  laser  diode.  The 


> 


laser  pulses  are  produced  at  the  radar  PRF  and  they  are  synchronized  to  occur  when 
the  acoustic  signal  in  the  AOD  consists  of  a  single  chirp  generated  by  the  SAW  device 
The  duration  of  each  laser  pulse  is  shorter  than  the  inverse  of  the  bandwidth  of  the 
chirp  signals  and  consequently  the  motion  of  the  travelling  acoustic  wave  does  not 
affect  the  operation  of  the  processor.  A  single  mode  AlGaAs  laser  was  used  in  the 
experiment,  at  15mW  peak  power  and  30  nanoseconds  pulse  width.  The  chirped  signal 
forms  a  phase  grating  in  the  AOD  with  linearly  increasing  spatial  frequency.  This  is 
shown  in  Fig.  5.  The  light  diffracted  near  the  piezoelectric  transducer  of  the  AOD 
in  Fig.  5  is  diffracted  at  a  larger  angle  than  the  light  diffracted  from  the  top  of 
the  device.  Consequently,  the  diffracted  light  comes  to  a  focus.  The  focal  plane  is 
located  11.5  cm  away  from  the  AOD  in  the  experiment  and  at  this  point  the  diffracted 
light  is  well  separated  from  the  zero  beam  and  it  can  be  detected.  The  self-focusing 
of  the  diffracted  light  accomplishes  the  focusing  of  the  SAR  image  in  the  range  (t) 
direction.  The  optical  field  at  the  plane  where  the  diffracted  light  comes  to  focus 
is  given  by 

S^(x,t)  «  J  S  (x’/v,t)  exp  [-jb^  (x’+x)^  /v2]  dx’ 

=  Sine  [  (x-2R  v/c)  T  b.  /ttv  ]  exp  [ j co  v  ( t-t  )  /cR  +  ltd  x],  (4) 

o  ol  o  a  v  o  o  o 


where  x'  and  x  are  the  spatial  co-ordinates  in  the  plane  of  the  AOD  and  the  detector 
respectively,  v  is  the  acoustic  velocity.  The  range  focused  signal  is  phase 
modulated  in  t  by  the  azimuth  chirp.  This  phase  modulation  can  only  be  detected 
interf erometrically .  Interferometric  detection  is  accomplished  by  adding  a  60  MHz 
sinusoid  to  the  radar  signal  and  applying  the  sum  to  the  AOD.  The  sinusoid  will 
cause  a  plane  wave  to  be  diffracted  in  the  same  direction  as  the  focusing  signal 
beam.  The  light  intensity  at  plane  x  is  given  by 


2 

^(x.t)  =  |si  (x , t )  +  exp  (jed  x)  I 

1  +  sine2  [(x-2R  v/c)  X  b_ /ttv2]  +  2  sine  [  (x-2R„v/c)  x  b./-v~] 

O  0  i  J  i 

X  COS  [(d  V  2  (t-t  )2  /R  c],  (5) 

o  a  o  o 


The  SAR  image  is  focused  in  the  azimuth  (t)  direction  by  correlating  the  time 

variation  of  the  above  function  with  the  function  cos  (cd^  v  t^/R  ).  This  is 

o  a  o 

accomplished  with  a  time  integrating  optical  correlator,  a  simple  diagram  of  which  is 
shown  in  Fig.  6.  A  transparency  with  intensity  transmittance  ^(y)  *  1  +  cos  (b^y^) 

is  illuminated  with  light  whose  intensity  is  temporally  modulated  by  I^(x,t).  y  is 

the  spatial  co-ordinate  perpendicular  to  x.  The  light  is  detected  after  the 
transparency  on  a  2-D  detector.  For  the  purpose  of  this  discussion  we  can  think  of 
the  detector  as  photographic  film.  If  the  film  travels  with  velocity  v'  in  the  y 


direction  while  it  is  exposed  to  light,  the  detected  signal  on  the  film. is 

I,(x,y)  ■  /  I.  (x,t)  I,  (y4v't)  dt.  We  choose  the  parameter  b,  ■  oi  v  /R  v'2  and 
-5  i  fc  z  o  a  o 

the  signal  accumulated  on  the  detector  becomes  the  autocorrelation  of  the  linear  FM 

waveform  in  Eq.  (5). 

I3  (x,y)  =  1  +  sine2  [(x-2Rq  v/c)  tq  b^/irv2] 

+  2  sine  [(x-2R  v/c)  t  b, /ttv2]  sine  [  (y-t  v')  b  Y/tt]  (6) 

ooi  o  2 


where  Y  is  the  total  distance  the  moving  detector  travels. 

The  first  term  in  Eq.  (6)  is  a  constant  bias,  the  second  is  a  signal  dependent 
bias  term  and  the  third  term  is  the  focused  SAR  image.  The  image  can  be  separated 
from  the  bias  terms  using  one  of  several  possible  methods.  For  instance  the  bias 
alone  can  be  calculated  separately  on  a  different  detector  and  then  substracted 
from  I^Cx.y)  to  produce  an  unbiased  image.  Alternatively  the  signal  term  in  Eq.  (6) 
can  be  formed  on  a  sufficiently  high  spatial  frequency  in  either  x  or  y  and 
separated  from  the  bias  terms  by  electronic  filtering. 

The  detector  used  in  the  processor  is  a  2-D  CCD.  The  photogenerated  charge 
stored  in  a  CCD  can  be  made  to  travel  along  one  of  the  spatial  dimensions  by 
applying  pulse  waveforms  to  the  device.  The  charge  travels  in  discrete  steps  but 
this  is  actually  convenient  for  this  system  since  t  is  a  discrete  time  variable  as 
well.  The  motion  of  the  photogenerated  charge  is  entirely  equivalent  to  the  physical 
movement  of  the  detector  in  Fig.  6.  Therefore  the  CCD  performs  the  correlation  in 
the  azimuth  direction  by  time  integration.  The  appropriate  value  for  the  velocity  v' 
is  set  by  transferring  the  change  pattern  on  the  CCD  by  one  pixel  at  the  radar  PRF. 

The  CCD  used  in  our  experiments  has  512  pixels  in  the  direction  of  charge 
motion  and  320  in  the  other  dimension.  The  entire  SAR  processor,  shown  in  Fig.  7,  is 
constructed  by  cascading  the  space  integrating  range  processor  of  Fig.  5  and  the  time 
integrating  processor  of  Fig.  6.  A  photograph  of  the  experimental  set-up  is  shown 
in  Fig.  8.  The  mask  used  as  the  azimuth  reference  is  shown  in  Fig.  9.  This  mask  was 
generated  with  a  computer  controlled  film  recorder.  Notice  that  the  chirping  rate  in 
the  vertical  direction  is  dependent  on  the  horizontal  (range)  position.  This  is 
necessary  for  compensation  of  the  range  azimuth  coupling  (the  focal  length  along 
the  azimuth  direction  of  the  zone  plate  in  Fig.  4  depends  on  the  range  of  each 
target).  A  photograph  of  the  focused  SAR  image  produced  in  real  time  with  the 
processor  in  Fig.  8  from  the  signal  displayed  in  Fig.  4  is  shown  in  Figure  10a. 

The  image  obtained  is  that  of  the  single  point  scatterer.  The  dimensions  of 
this  focused  spot  along  the  horizontal  (r&nge)  and  vertical  (azimuth)  directions 
determine  the  resolution  of  the  SAR  processor.  The  width  of  the  spot  in  the 
horizontal  direction  is  expected  to  be  equal  to  the  acoustic  velocity  divided  by  the 
bandwidth  of  the  chirp  (20  MHz)  produced  by  the  SAW  device  [3].  The  width  of  the 
spot  in  Fig.  10a  is  broader  than  expected  by  a  factor  of  two.  We  attribute  this  to 
apodlzatlon  due  to  the  gausslan  shape  of  the  wavefront  Illuminating  the  AOD  and 
possibly  abberatlons  in  the  lenses.  The  number  of  pixels  that  the  spot  is  expected 
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to  occupy  in  the  vertical  (azimuth  direction)  is  equal  to  the  total  number  of  pixels 
of  the  CCD  in  the  vertical  direction  (512)  divided  by  40  (the  time-bandwidth  product 
of  the  azimuth  chirp).  The  measured  spot  width  is  approximately  equal  to  what  is 
expected  in  the  azimuth  direction.  The  Isometric  display  of  the  focused  image 
(shown  in  Fig.  10b)  reveals  the  structure  of  the  bias  on  which  the  image  is  formed. 
There  is  a  flat  plateau  due  to  the  collimated  reference  wave  and  for  this  case  of  a 
single  point  target  an  additional  ridge  along  the  azimuth  direction  (the  second  term 
in  Eq.  3).  The  sidelobe  structure  of  the  focused  spot  can  be  observed  in  the  cross 
sectional  scan  (along  the  azimuth  direction),  shown  in  Figure  10c.  In  Figs.  10b  and 
10c  the  compression  in  the  azimuth  direction  is  80  (2  times  more  than  in  Fig.  10a). 

The  digital  simulation  of  the  autocorrelation  of  a  linear  FM  waveform  with  a  time- 
bandwidth  product  equal  to  80  (displayed  in  Fig.  40d)  shows  that  the  level  and  the 
structure  of  the  optically  produced  sidelobes  are  remarkably  accurate.  The  sidelobe 
level  can  be  reduced  further  by  apodization  and/or  increasing  the  pulse  compression 
ratio  in  the  azimuth  direction.  The  value  80  for  azimuth  pulse  compression  we  used 
in  the  experiment  is  typical  for  aircraft  radars. 

ADAPTAELE  ARCHITECTURE 

In  some  applications  the  parameters  of  the  radar  such  as  the  velocity  and  the 
altitude  of  the  aircraft  and  the  direction  in  which  the  antenna  is  pointed  (side 
looking,  forward  looking,  spotlighting)  change  dynamically.  Consequently,  the  real 
time  SAR  imager  must  have  the  capability  of  adapting  rapidly  to  such  changes  in  order 
to  continuously  produce  a  well  focused  image.  The  impulse  response  of  the  processor 
described  in  the  previous  section  is  determined  in  part  by  the  transmittance  of  the 
mask  used.  A  real  time  spatial  light  modulator  can  be  used  in  place  of  the  mask. 

Since  the  transmittance  of  such  a  device  can  be  changed  dynamically,  the  impulse 
response  of  the  processor  can  be  adapted  to  match  the  radar  parameters.  An  alternate 
and  preferable  approach  is  to  use  a  second  AOD  in  place  of  the  fixed  mask.  The  mask, 
as  we  see  in  Fig.  9,  introduces  a  2-D  spatial  light  modulation.  The  modulation  along 
one  of  the  dimensions  however  is  simply  scaling.  Therefore  the  AOD,  which  is  a  1-D 
spatial  light  modulator,  can  be  used  to  introduce  the  basic  chirp  modulation  in  the 
y  direction.  An  optical  system  that  images  the  AOD  along  the  y  direction  with 
different  magnification  for  each  x  position  produces  a  2-D  spatial  light  distribution, 
which  when  detected  interferometrically ,  results  in  a  pattern  identical  to  the  one 
shown  in  Fig.  9.  This  is  shown  schematically  in  Fig.  11.  A  cylindrical  lens  is 

used  to  focus  the  light  in  the  AOD  and  recollimates  the  diffracted  light.  The 

lens  is  a  cylinder  that  has  focusing  power  in  the  vertical  direction.  Let  us 

assume  that  the  AOD  is  modulated  by  a  single  frequency.  The  diffracted  light  in  this 
case  is  a  plane  wave  and  it  will  be  focused  to  a  horizontal  line  by  the  cylindrical 
lens  Ly  is  rotated  around  a  pivot  coinciding  with  the  y  axis  by  an  angle  0. 

The  focused  line  lies  in  the  (x,z)  plane  and  it  is  at  an  angle  0  with  respect  to  the 
x  axis.  The  light  propagates  as  a  diverging  tilted  cylindrical  wave  after  the  focal 
plane  of  L^.  If  we  observe  interferometrically  this  light  distribution  in  a  plane  ?2 

perpendicular  to  the  z  axis,  we  will  see  the  pattern  shown  in  Fig.  4.  If  the  inter¬ 
ference  of  the  output  of  the  systems  in  Figs.  5  and  11  is  formed  on  the  CCD,  the 
detected  signal  will  be  exactly  the  same  as  the  signal  in  Eq.  5.  It  is  interesting 
to  note  that  the  tilted  cylindrical  lens  alone  is  sufficient  to  produce  the  desired 
interference  pattern.  The  AOD  however  provides  added  flexibility  in  several  respects. 
The  chirping  rate  and  the  starting  frequency  of  the  azimuth  reference  can  be 


dynamically  controlled  by  modulating  the  AOD  with  the  appropriate  temporal  waveform. 
Furthermore,  the  azimuth  correlation  can  be  formed  by  shifting  the  signal  in  the 
AOD  instead  of  transferring  the  charges  in  the  CCD.  This  capability  is  useful  for 
producing  multiple  looks,  i.e.,  forming  independently  several  images  of  the  same 
area  and  adding  them  to  reduce  speckle.  If  the  charge  pattern  in  the  CCD  is 
travelling  during  exposure,  the  focused  images  are  continuously  read-out  and  there¬ 
fore  it  is  not  possible  to  add  on  the  detector  multiple  looks.  Finally  the  use  of 
the  AOD  allows  us  to  build  the  architecture  as  a  cascade  of  the  systems  in  Figs.  5 
and  11,  rather  than  a  two-arm  interferometer.  This  is  an  important  consideration 
when  we  consider  the  mechanical  stability  of  the  system.  The  cascaded  architecture 
is  shown  in  Fig.  12.  In  this  architecture  all  beams  go  through  the  same  optics. 

In  order  to  record  a  chirp  signal  on  the  detector  we  must  interfere  two  waves  with 
different  degrees  of  collimation.  The  second  AOD  in  Fig.  12  is  modulated  by  the  sum 
of  a  single  sinusoid  and  a  chirp  signal.  The  light  diffracted  by  this  AOD  consists 
of  a  plane  wave  and  a  diverging  spherical  wave  and  thus  the  interference  of  these 
two  diffracted  waves  after  lens  produces  fringes  with  linearly  increasing  spatial 
frequency  in  the  y  direction,  as  desired. 
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Figure  3.  Quad-  ahase  modulation  simulating 
’ '  -  ..h  phase  history. 


Figure  4.  Simulated  SAR  impulse  response  produced  in 
real  time  by  the  system  in  Figure  2. 


Figure  5.  Acousto-optic  processor  that  performs  the 
focusing  of  the  SAR  image  in  range. 
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Figure  6.  Time  integrating  optical  proceaaor  that  performs 
the  focusing  of  the  SAR  image  in  aximuth. 


Figure  9.  Computer  generated  mask  used  in  the  processor 
of  Figure  8. 


Figure  10«.  Focused  SAR  image  of  the  single  point  scatterer 

corresponding  to  the  zone  plate  in  Figure  4,  produced 
ia  real-time  with  the  processor  of  Figure  8. 
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Figure  10b.  Isometric  display  of  the  focused  point  scatterer  demonstrating 
structure . 

10c.  Cross  sectional  scan  in  azimuth  shoving  the  sidelobe  structure 
compression) 

lOd.  Digital  simulation  of  the  signal  shovln  in  Figure  10c. 
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Figure  11.  Acousto-optic  system  with  which  the  azimuth  reference  is  interfer- 
ometrically  introduced. 


•court  (Mrric 
areas 


“  As* 

10^1 

msm  I 


PfTO 

SIMM. 


Figure  12. 


Optical  architecture  for  SAK  imaging  with  the  azimuth  reference 
introduced  through  a  second  ADD  In  place  of  the  mask  used  in  Figure  7. 
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REAL-TIME  SYNTHETIC  APERTURE  RADAR  IMAGE 
FORMATION  UTILIZING  AN  ELECTROOPTIC  SPATIAL 
LIGHT  MODULATOR 
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and 
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Jet  Propulsion  Laboratory 
California  Institute  of  Technology 
Pasadena,  CA  91109 


In  the  process  of  synthetic  aperture  radar  data  acquisition, 
each  point  target  in  the  terrain  to  be  imaged  is  illuminated  by 
numerous  consecutive  radar  pulses  during  an  azimuthal  terrain 
traverse  by  the  radar  platform  (typically  aircraft  or  space 
vehicles)  .  Thus  a  given  point  target  contributes  a  return  signal 
that  mimics  the  radar  pulse  shape,  at  a  time  delay  that  depends 
on  the  point  target  range  at  each  radar  platform  azimuth,  and 
at  a  center  frequency  Doppler-shifted  by  the  instantaneous 
relative  velocity  of  the  platform  with  respect  to  the  target. 
Recovery  of  the  point  target  location  and  scattering  strength 
from  such  a  two-dimensional  (range  and  azimuth)  signal  record 
is  exceedingly  complex,  since  the  return  signal  amplitude  at 
at  any  given  time  is  comprised  of  scattered  components  from 
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all  points  in  the  illuminated  terrain. 

The  image  formation  algorithm  from  synthetic  aperture 
radar  data  involves  the  two-dimensional  deconvolution  of  the 
recorded  phase  and  amplitude  history  of  the  radar  return  from 
the  radar  pulse  shape  in  the  range  dimension,  and  from  the 
Doppler-induced  phase  modulation  in  the  azimuth  dimension. 

For  ease  in  digital  processing  using  fundamentally  sequential 
architecture  machines,  the  two-dimensional  deconvolution  is 
typically  replaced  by  two  one-dimensional  deconvolutions  in 
the  azimuth  and  range  dimensions.  The  reference  function 
utilized  in  the  azimuth  deconvolution  is  a  function  of  range, 
and  hence  is  inherently  space-variant.  Since  the  reference 
functions  employed  typically  require  of  order  1000  sampled 
elements  to  maintain  high  resolution  in  the  reconstructed 
image,  such  computations  are  quite  cumbersome  and  often 
prohibitive . 

Optical  reconstruction  of  synthetic  aperture  radar  images 
can  be  readily  accomplished  by  recording  the  return  signal 
phase  history  on  an  image  storage  device  such  as  photographic 
film.  The  two-dimensional  amplitude  that  results  from  a  given 
point  target  forms  a  zone  plate  that  is  self-focusing  in  both 
dimensions  when  coherently  illuminated.  Since  the  azimuth  and 
range  focal  lengths  are  functions  of  distinct  parameters,  they 
must  in  general  be  brought  into  coincident  focus  by  cylindrical 
lenses  oriented  to  act  separately  on  the  range  and  azimuth 
dimensions . 

A  number  of  advantages  accrue  to  such  optical  SAR  processors 


due  to  the  inherent  parallelism  of  the  optical  processing 
algorithm.  The  inherent  speed  advantage  allows  real-time 
systems  to  be  envisioned  that  are  small,  lightweight,  and  low- 
power.  Although  digital  processors  are  usually  characterized 
by  superior  computational  accuracy,  they  are  typically  quite 
large,  expensive,  and  power  consumptive.  Many  radar  applications 
such  as  quick-look  or  search  modes  require  or  would  be  better 
served  by  on— board  processing.  These  SAR  systems  will  require 
small,  lightweight  processors,  especially  for  use  in  conjunction 
with  spacecraft-borne  radars. 

Present  implementations  of  coherent  optical  SAR  processors 
share  a  common  limitation,  in  that  real-time  operation  is 
obviated  by  the  use  of  photographic  film  as  an  incoherent- to- 
coherent  image  transducer.  One  possible  approach  to  high  speed 
optical  synthetic  aperture  radar  signal  reconstruction  involves 
the  utilization  of  two-dimensional  real  time  spatial  light 
modulators  as  recycleable  replacements  for  photographic  film  in 
the  input  transducer  plane  of  modified  SAR  coherent  optical 
processors,  as  shown  in  Figure  1.  Leading  candidate  spatial 
light  modulators  include  the  Pockels  Readout  Optical  Modulator 
(PROM),  CCD-addressed  liquid  crystal  light  valves,  and  CCD- 
addressed  membrane  light  modulators.  The  fundamental  physical 
limitations  of  such  real  time  devices  that  affect  SAR  processor 
performance  characteristics  are  under  investigation.  Current 
research  on  the  PROM  is  focused  on  the  effects  of  device 
operational  mode,  device  constitutive  parameters,  electrooptic 
crystal  orientation,  writing  wavelength,  frame  rate/data  overwrite/ 
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presuming,  erasure  completeness  and  image  retention  on  the 
overall  quality  of  SAR  image  formation. 

Several  modes  of  addressing  such  tv.’o-dimens ional  spatial 
light  modulators  are  under  active  consideration,  including  a 
z-axis  modulated  x-y  laser  scanner,  an  image-intensified  CRT 
display,  and  an  acoustooptic  delay  line.  A  promising  new  device 
the  Xerox  Linear  Array  Total  Internal  Reflection  spatial  light 
modulator,  is  under  development  and  evaluation  for  the  page 
composition  function. 

The  recent  development  of  such  novel  electrooptic  spatial 
ight  modulators  (both  one-and  two-dimensional)  will  allow  real¬ 
time  coherent  optical  SAR  processors  to  be  implemented.  It  is 
hoped  that  the  availability  of  such  processors  will  initiate 
new  mission  applications  not  presently  envisioned  due  to  current 
digital  processor  limitations. 
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